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as alternative du principe de focalisation par réflexion et du principe 


e tube Dent occuper Pune. ou d'autre de due positions nu node 
da du couvercle de la boîte; corps et anticathode solidaires du fond. Le réglage 

d à ch ue cas par rotation du spectrographe autour d’un axe vertical. passant DATA 

os et au au de focalisation. La plaque PHotogrephidues se ae 


qui dé déjà ere de à. 18 À 2! les nt 
sont très” satisfaisants. Je vais. donc décrire ce. 
nouveau montage. d 


FA 
Pa 


| )Principe. — Le So ot utilise comme 
ue pe ik mous. “Nous sentions analyseur une lame cristalline courbée suivant un 
spectrographe à vide pour étendre cylindre; cette lame travaille selon deux, principes 
TS. les Mu ne de grandes _ de focalisation CONNUS : 

Pour l'exploration des petits. Anbles de Bragg 
. (donc, conformément à la formule de Bragg : : 
nÀ = 2 d sin +, pour l'exploration des petites lon- 
puisse ne dde d’un  gueurs d’ onde), le spectrographe réalise le principe 
équences très étendu : depuis les plus de focalisation que j'ai développé [x ] pour les fais- 
de ee jusqu’à la  ceaux réfléchis sélectivement sur les plans réticu- 
laires internes d’un cristal; la lame travaille par 
transmission. La source de rayons X doit être placée 
extérieurement au spectrographe proprement dit. 
J’appellerai ce cas le cas T (transmission). 

Pour l'exploration. des grands angles de Bragg, 
c’est-à-dire des rayons X mous (ainsi que des rayons 
durs ou moyens dans des réflexions d'ordres supé- 
ricurs), le spectrographe réalise le principe de foca-: 
lisation indiqué par Johann [2]. La réflexion sélective 
a lieu alôrs sur la face cylindrique concave de la 
cours l'année no due 1943- a J’en ai lame courbée. La source de rayons X est placée, 
l'épreuve dans une zone de longueurs d'onde par rapport à cette surface, du même côté que 
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l'axe même du ‘cylindre de courbure. a appellerai 
ce cas le cas R (réflexion). 

J'ai rappelé ces deux principes de nord X, 
désormais classiques, dans les deux schémas ci- 
contre (fig. 1). La focalisation qui peut être consi- 
dérée comme sensiblement parfaite pour les rayons 
centraux Correspondant aux ouvertures pratique- 
ment possibles pour le cristal courbé, produit un 
élargissement des raies spectrales dont rendent 


compte les rayons extrêmes que j'ai tracés sur 
chaque schéma. 


CasT 
RISE. 


Le figures montrent bien que la disposition 


respective de la lame cristalline cylindrique et du 


châssis photographique destiné à l’enregistrement 
des spectres, ne doit pas différer essentiellement dans 
les cas T et R. Si bien qu'un même assemblage, 
judicieusement prévu, devait pouvoir constituer, dans 
les deux cas, le spectrographe proprement dit. 
Nous avons décrit ailleurs [3], [4] le modèle de 
spectrographe que nous avons trouvé le plus pra- 
tique, entre divers autres, pour travailler à l’air en 
utilisant la focalisation par les cristaux courbés. 
Le cylindre de focalisation y est matérialisé par un 
cylindre métallique; le châssis photographique. est 


mobile sur ce cylindre auquel la plaque reste tan-. 


gente; le porte-cristal s'appuie contre sa partie 
intérieure, par sa propre base qui est tournée au 
rayon exact. La figure 2, extraite d’un précédent 
article [4 reproduit le schéma de l'appareil, schéma 
adopté, ici encore, à peu de choses près. 

L'étude des rayons X mous devient pratiquement 
impossible dans. l’air au delà de 2 ou 3 À. L’utili- 
sation du spectrographe pour les plus grandes 
longueurs d'onde exigeait donc qu’il soit enclos dans 
une enceinte vidée. Il est indiqué d’y comprendre 
aussi le tube à rayons X. IL y a intérêt, spécialement 
pour les conditions de technique du vide, à donner, 
à cette enceinte, des dimensions aussi petites que 
possible. 
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cas R se ferait par un changement de position du 


Ces Gel eee points 
s'agissait de trouver, pour notre montage, u ine 
sition avantageuse du tube et. du spectrographe 
Nous avons décidé que le passage du cas T au 


tube dans l’enceinte, tandis que le réglage du spec 
trographe par rapport au tube se ferait dans chaque 
cas, ainsi que dans nos montages fonctionnant à 
l'air, par rotation du spectrographe autour d’un axe 
vertical qui coïncide avec l’axe de tangence du 
cristal et du cylindre de focalisation. Nous nous 
sommes ainsi arrêtés au principe suivant (fig. 4) : 

L'’enceinte vidée est un grand cylindre d'acier 
dont la base horizontale Best fixe sur un bâti 
massif auquel sont assemblées les pompes à diffu- 
sion et e canalisations de HORIDASS: Le « couvercle » » 


#1 


Rotation du 
‘spectrographe 4 
sur sa table, : 


un Fixation da 
spectrographe. 


Spectrographe à rayons K'atilisant des lames cristallines courbées, 
par transmission, (Schéme du principe de la construction). 

B: Coupe verticale re 

B—Pb. D—Porte cris- 
tal. E—lame cristal- 


À: Projection horizontale 


B— Ecran de Pb. D— Porte cristal smovible. 
G—Graduation indiquant la position du 


spectrographe. H—Visour à pointe pour line. I—Châssis 
les réglages. I— Châssis photographique photographique 
mobile, mobile amovible. 


J— — Plaque photo. 
K—Table ‘portant 
l'appareil, 


Fig. 2 


horizontal C peut être soulevé, à l’aide d’un système 
à roue dentée et crémaillère, parallèlement au 
plateau B, soit seul, soit solidairement avec la paroi 
cylindrique P. Le spéctrographe même repose 
sur B; il peut tourner, autour de l’axe À, de plus 
de 459 de part et d’autré de sa position moyenne, 
lorsqu'il fonctionne par transmission. Le tube à. 
rayons X prend place soit en T, soit en R, suivant 
que le cristal est utilisé par transmission ou par 
réflexion. La direction TR passe par À; elle fait un 
angle de 30° avec l’axe diamétral du spectrographe | 
lorsque celui-ci occupe sa position moyenne. Les … 
distances AT et AR sont égales (fig. a à 
Certaines des conditions précisées n'étaient pas # 
indispensables; elles ont été choisies pour des raisons » | 
de commodité. Par exemple, le choix des positions T è 
et R pour le tube permet, quelle que soit celle de } 
ñ 
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Exemple de spectres obtenus par transmission, 
Mica (201). Région spectrale : 0,8 à 1,5 À. 


(1) 
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(3) 
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%y, le XX, 
nm 
Re X2 Ca.K 10€ordre 
Exemple de spectres oblenus par réflexion, avec le même mica. — 
Anticathode : Al, Temps de pose croissant de (1) à (3). Région spectrale : 16,8 À. 


Photo 1. 


Photo 3, 


Photo 2. 


“4 Le pr b: a pur assure, ons le corps du 

tube, un pompage plus rapide et plus parfait que 
- dans l’ensemble de l'enceinte. Par ailleurs, grâce 
_ à des systèmes de joints et fixations interchan- 
ph le tube S ‘adapte, DA nent en T ou 


arties î + “Corps: auquel Rue eaihode. 
_est solidaire du plateau B. L'isolant, qui porte la 
_cathode, repose sur le plateau supérieur C. Lorsque 
la grande boîte BCP est fermée, le tube se 
uve constitué : l'extrémité supérieure du corps 
: vient s'appuyer contre la face interne du plateau C 
par. l’intermédiaire d’un joint de caoutchouc et la 
cathode prend sa place en face de l’anticathode. 
IL faut évidemment qu'ait été réalisé, par construc- 
: tion et assemblage, un centrage. parfait des diffé- 


. Lorsque le jube occupe la one T. les deux 
ouvertures libres (respectivement dans B et dans C) 


mesure du vide; ces ajutages sont placés dans les 
ouvertures T lorsque le tube est monté en R. 

Afin que le corps du tube puisse occuper cette 
position R, il est nécessaire de forer dans le plateau 
_ du spectrographe lui-même une ouverture qui le 
<4 
_ laisse passer librement et assure la possibilité d’une 
rotation du spectrographe autour de.l’axe À dans 
un domaine angulaire assez étendu. J’ai fait pra- 
“à tiquer une ouverture qui n’entame que le plateau 
- et respecte, presque intégralement, le rebord cylin- 
 drique du spectrographe; ainsi, la rigidité reste 

_ parfaite. Étant donné le diamètre du corps du tube, 


1 


1 , . . . pale > DT # 
“ ration une limite inférieure égale ici à environ /o0. 


è 


4 _ ceci impose aux angles de Bragg susceptibles d'explo- 


Cette limitation n’est pas grave, car les conditions 
de réglage et de bonne focalisation, dans la méthode 
de Johann, rendent peu avantageux l'usage du 
… spectrographe pour des angles plus petits. Les rayon- 
… nements correspondant à ceux-ci s’analysent plus 
‘à favorablement par transmission. 

D  L'échancrure pratiquée dans le spectrographe 
3 s'étend, d'autre part, jusqu'au voisinage du rayon 
médian. Elle permet d'atteindre le plus grand angle 
de Bragg pour lequel le faisceau réfléchi par la face 
concave du cristal puisse frapper la piaque photo- 
graphique sans être intercepté par le corps même 
: du tube. Cet angle est égal ici à 81° 30’. 


. Caractéristiques du montage réalisé. Ee 
” diamètre du spectrographe est de 25 cm. Les porte- 
> cristaux sont des « lentilles cylindriques convexes » 


ns UNIVERSEL » POUR LES 


‘déplace le spectrographe, 


rentes pieces et que ce centrage se conserve après 


RAYONS ù Sr it .91 


‘48 5o cm de diamètre. Cette. dimension correspond 


aux possibilités extrêmes du tour-dont nous dispo- 
sions ét tombait bien dans l’ordre de grandeur que 
nous avions choisi. Nous avons dû, en conséquence, 
donner à la boîte vidée où se fixe le tube et où se 
un diamètre intérieur 
de 5r0o mm. Les schémas (fig. 3, 4 et 5) et les photo- 
graphies (au verso de la planche) montrent la dispo- 
sition des différentes pièces et l'aspect du montage 
réalisé. 


La grande boîte vidée est en acier ordinaire: Le 
spectrographe et le tube sont en acier inoxydable. 
Voici encore quelques détails complémentaires : 


Dispositif de pompage. — Le vide est assuré par 
un groupe de deux pompes métalliques, la première 
à diffusion de mercure, la seconde à diffusion d’huile. 
Les caractéristiques indiquées par le constructeur 
(C. GR.) sont les suivantes : vitesse de pom- 
page, 351; vide limite, < 10 *baryes. La jauge 
ionisation (C. G.R.), montée sur l’enceinte, indique 
des courants résiduels de Flordre de 0,3 A qui 
correspondent à peu près à 10 ° mm Hg. 

Le vide primaire comporte une pompe rotative 
à palettes, à deux étages, qui peut être mise en 
liaison, à travers un réservoir vidé auxiliaire, avec 
un relai connecté à une résistance calibrée fonc- 
tionnant en jauge thermique de contrôle; de cette 
manière, la pompe à palettes peut n'entrer en 
service que lorsque la pression dans le réservoir- 
tampon s'élève au-dessus d’une limite fixée. 


Entre b pompe à diffusion d’ huile et les canali- 
sations allant à l’enceinte et au tube à rayons X, 
un gros robinet métallique à trois voies est inter- 
posé (N. B.). La rotation d’un axe extérieur entraîne, 


à travers un petit rodage conique, le mouvement d’une 


crémaillère intérieure. Cette crémaillère, qui se 


À 


Fig. 4. — Montage pour la transmission (cas T). 


trouve donc placée sous vide, commande deux 
petits pistons qui ferment alternativement, par 
compression de joints de caoutchouc contre des 


épaulements convenables, l’une ou l’autre extrémité 
d’une tubulure de liaison dont le diamètre intérieur 


est ici de 5o mm. Lorsque la crémaillère est dans la 
position de « fermeture », le piston correspondant 
assure l’obturation de la sortie de la pompe; simul- 
tanément, l'enceinte spectrographique se trouve en 
liaison directe avec la pompe rotative, par une petite 
canalisation latérale. Lorsque la er‘maillère est dans 
la position « d'ouverture », la pompe à diffusion 
d'huile est mise en communication avec l'enceinte, 
tandis que, par là même, la connexion de cette 
dernière avec la pompe à palettes est interrompue. 

Ce dispositif présente des avantages intéressants : 


il joue le rôle que jouerait un gros robinet ordinaire 


à clef et boisseau coniques rodés et graissés, sans en 
présenter les inconvénients graves connus : limi- 
tation de la vitesse de pompage, tension de vapeur 
des graisses, risque de « gripage » et de fuite, etc. 
De plus, la communication auxiliaire directe entre 
le spectrographe et la pompe primaire permet de 
laisser les pompes à diffusion en fonctionnement 
continu, même si l’on doit ouvrir le spectrographe 
pour y effectuer réglages ou changement de plaque 
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. deux cylindres d’acier inoxydables entre lesquels 


Dhotographique uffit de 
dans l'enceinte avant de 
la pompe à diffusion d'huile. 
conduit à un gain de temps cond tale 
Les canalisations de vide sont toutes montées 
l'intermédiaire de joints de caoutchouc. L’ lasticit 
est obtenue, s’il en est besoin, par des tubulur 
métalliques souples : tubes « tombac » et surt L 
membranes thermostatiques. 


Il est important nes ne . laisser la Pts 


n’est donc ouvert que M la pompe à diffu 10 
d'huile est déjà assez chaude pour que son. je 
empêche la distillation du mercure. au ARE je ilr 


(1 


# 


de mercure sur PRaticatpode qe cool 
tement toute autre émission X caractéristique, d 
grande longueur d’onde. Pour la même raison, 1 
vide, dans de tube proprement dit, doit être aus 
poussé que possible, afin d'éviter au maximum 1 
pulvérisation cathodique et la formation de dépôt 
par calcination des graisses ou huiles résiduelles sur 
l’anticathode. D'où l'usage d’une canalisation spé 
ciale de pompage, sur laquelle, dans certains cas, 
il est bon de prévoir une zone de condensation su 
air liquide. 


Le tube à rayons X. — Le corps est constitué pars 
est établi un courant d’eau pour son refroidissement. ” 


À toutes fins utiles, la fenêtre de sortie du rayon- 
nement a été forée assez large : 13,5 sur 5 mm. 


a “A sis est bréfé- 
‘toujours masquer l'orifice de sortie du 


A. 


très mince, puisque la différence de pression 
entre ses deux faces est négligeable. Toutefois, 
A est bon que sa densité superficielle soit suffisante 
pour qu'elle arrête les électrons diffusés, réfléchis par 
icathode, ou de toute autre origine, qui pour- 
ent, sans cela, franchir l’orifice et produire des 
s pernicieux à l’extérieur du tube. 
s très mous, il faut utiliser des feuilles d’alu- 
um battu ou, mieux encore, de la _cellophane 
mince de couleur rouge pour masquer LEURS, 
lumineuse de la cathode. ss 
‘isolant est un cylindre de stéatite de à 
j chaque extrémité, par une joue 
plane. Les joints de caoutchouc, 


mitoyennes, viennent s’écraser 
contre ces joues 1 5 
La pièce cathodique comporte 
eee ---Q un raccord métallique souple qui 
“permet d’en réaliser le centrage. 
Celui-ci, une fois effectué le 
changement de cathode, se fait 
en démontant la pièce au-dessus 
de ce raccord, sans plus y tou- 


cher. Il est fait usage de 
cathodes chaudes : filaments 
\ « reCtilignes » ou cathodes à 


oxydes. 


passé 1 kW, il n’a pas été utile 
de refroidir la tête de cathode. 

La cathode est logée au sein 
d'une cupule de concentration. 


gène et, dans le cas R, assez 
hard pour que le spectre soit 
assez étendu. On agit en conséquence sur la posi- 
ion de la cupule et de la cathode. 
si Nous avons établi plusieurs 
_cathodes : 
L'un comporte une longue base refroidie, en 
_ cuivre rouge, dont la tête est filetée. Sur ce filet 
se peuvent se visser l’une quelconque des nombreuses 
#1 têtes. d’anticathodes interchangeables dont nous 
poor au laboratoire. Nous avons ainsi immé- 
 diatement un jeu d’éléments anticathodiques à l’état 
massif. Les autres éléments doivent être fixés sur 
une tête de cuivre ou d'aluminium, par exemple. 
. Une autre anticathode, en cuivre massif, comporte 
trois faces. L'une ou l’autre de ces faces peut être 
amenée dans l’axe du tube, vis-à-vis de la tête de 


modèles d’anti- 


24 Re 


x 


nement par une feuille mince qui réalise 
à dépendance complète du tube par rapport à. 
4 nsemble de l'enceinte vidée. Cette feuille peut 


Pour les. 


ortés par les parties étalli ues . 
ge Ê à ti 7 deviennent, elles-mêmes, trop absorbantes. 


toute la hauteur 


Pour les puissances. 
consommées. qui n’ont pas dé-. 


supports par des bracelets de caoutchouc de diamètre 


La tache focale doit être homo- 


passage du rayonnement. Un même cristal peut 


pi Son constituer une anticathode. La 
manœuvre peut se faire sous vide, grâce à un rodage 
conique excentré, usiné à la partie intérieure de 120 
pièce (fig. 6). Ceci permet de photographier, sur 
un même cliché, sans effectuer de rentrée d’air, 
les spectres des éléments présents sur les trois faces. 


ss 


Spectrogr aphe. — Je n’ai plus que peu à dire sur 
le spectrographe lui-même. Il diffère cependant de 
ceux que nous employons à l’air par quelques détails 
qui tiennent à la nécessité de le placer sous vide. 

Ainsi, l'enregistrement des spectres de rayons X 
assez pénétrants se fait aisément à l’aide de châssis 
photographiques clos, étanches à la lumière ordi- 
naire; la surface sensible est protégée du jour par … 
une feuille de papier noir, que le rayonnement 
traverse sans affaiblissement notable. Si le châssis 
est destiné à être utilisé sous vide, il est préférable 
de remplacer le papier par de l’aluminium mince. 
Mais si, de plus, il s’agit d'observer des rayonnements | 
très mous, lés feuilles d’aluminium laminées les. 
plus minces, encore exemptes de perforations, 
Il faut 
offrir aux faisceaux réfléchis la surface sensible nue. 
Ceci oblige à rendre l'enceinte spectrographique 
étanche à la lumière et à ne manipuler qu'en chambre 
obscure. Il faut renoncer à faire usage de la jauge # 
à ionisation, car son filament est fortement lumi- 
neux. Pour obvier à ces inconvénients, nous construi- 
sons un châssis étanche au jour, müni de deux volets 
qui seront ouverts sous vide, au moment de prendre … 
le spectre, à l’aide d’une petite commande électro 
magnétique. Il sera possible d’enregistrer, soit. 
du spectre, soit seulement sa 
demi-hauteur, ce qui est commode, par exemple, 
pour superposer au spectre à mesurer des raies 
« de référence ». ma 

D'autre part, nous trouvons généralement excel. 
lente la fixation dés lames cristallines sur leurs 


convenable. Ce procédé n’est pas très recommandable 
pour les porte-cristaux qui doivent être placés sous 
vide. Il est alors préférable d'assurer le serrage du 
cristal par une contre-pièce métallique. J’emploie 
de minces lames de ressort en acier, élastiques et 
percées en leur centre d’une ouverture pour le 


servir indifféremment par transmission ou par 
réflexion s’il cst orienté par rapport à l’axe du 
cylindre de courbure pour offrir aux réflexions internes 
un système réticulaire intéressant. On sait que 
cette orientation du cristal avant courbure n’est 
pas nécessaire si l’on fait usage de la réflexion sur 
la surface de clivage. La méthode par réflexion 
semble exiger de meilleures qualités optiques de la 
surface cylindrique concave du cristal que la méthode 
par transmission, laquelle exige, par contre, une 
homogénéité plus parfaite du cristal dans toute son 
épaisseur. Le choix des lames doit donc être fait 


tante réticulaire d — 
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avec un soin d'autant plus rigoureux que leur emploi 


doit être double. Il est effectué, comme d'habitude, 
à la fois par des procédés optiques et aux rayons X. 
Il en est de même du contrôle de leur montage dans 
les porte-cristaux, contrôle que l’on termine en pre- 
nant uné série de spectres par lesquels il s’agit de 
délimiter la région d’une lame à peu près bonne 
qui donne des raies spectrales vraiment parfaites. 


Cette région-est isolée par des diaphragmes métal- 


liques mis en place une fois pour toutes. Une partie 
des contrôles peut se faire à l'air, en réglant le 
spectrographe devant un tube autonome, ce qui est 
plus rapide. 


Choix de la nature des cristaux. —— Sous le 


rayon de courbure de 25 cm, il est relativement facile 


de courber élastiquement une série de cristaux, 


comme le mica, le gypse, le quartz etc.; de courber 
plastiquement le sel gemme et d’autres. 


Mica. — Le mica a pour lui sa remarquable élas- 


ticité. Il offre, à la méthode de transmission, un 


certain nombre de systèmes réticulaires intéressants, 
parmi lesquels nous employons souvent les plans (201) 


et (100) qui, appartenant à la même zone, sont 


accessibles dans un même montage de la lame. Je 


rappelle leurs caractéristiques : pour (201), incli- 
maison sur la normale aux faces & = 40 56, COR 
2,65 À ; pour (100), & — 109 ro’, 
d = 2,61 À, pouvoir réflecteur très faible en deuxième 
ordre. L'un et l’autre donnent ici une dispersion 


. d’environ 20 ux/mm en premier ordre. 


_ Le mica n'offre, à la méthode par réflexion, que 
ses plans de clivage (001) dont le pouvoir réflecteur 
est élevé dans de nombreux ordres successifs; la 
constante réticulaire qui, ainsi que celles précé- 


- demment indiquées, varie quelque peu d’un genre 
de mica à l’autre, est à peine inférieure à 10 À (1). 


Il s'ensuit que la dispersion en premier ordre est 
faible, mais que les micas permettent l'étude des 
rayons X mous jusqu'à une vingtaine d’angstrôms. 
Pour aller plus loin, il faudrait faire usage de chlo- 
rites, pour lesquelles d,9, — 14 À. Mais il semble 
difficile de trouver des lames utilisables (?). Au delà, 
il n’y a plus de cristaux naturels, et il faut faire 
appel soit aux couches orientées d'acides gras, soit 
aux réseaux optiques. Ces rappels sont valables pour 
toute méthode de spectrographie X par cristaux 
dans ce domaine. 

Dans notre montage, comme je l’ai montré, les 
angles de Bragg, susceptibles d’observation, sont 
compris pour le cas T entre o et 750, et pour le cas R 
entre 4o° et 81° 30’. L'analyse de toutes les lon- 
gueurs d'onde qui relèvent de la méthode de trans- 


(:) On se reportera au traité de M. SEGBAHN, Spektroskopie 
der Rôüntgenstrahlen, 1931. 

(2) Je remercie M. d’Orsel d’avoir bien voulu me prêter 
quelques beaux échantillons de chlorites pour essais. 


peut être orientée avant courbure de manière à. 


mission est possible avec n'importe quel cristal. 

La limitation du domaine d’angles détermine le : 
domaine correspondant de longueur d'onde obser- 
vable avec notre appareil, grâce à une réflexion. 
cristalline partieulière. Pour les plans (001) du 
mica, les réflexions d'ordres successifs rendent obser- 
vables des gammes discontinues de longueurs 
d'onde qui ne se recouvrent pas. Le tableau ci-contre 
précise ce fait. L'emploi du mica ne suffit donc pas 
pour l’üniversalité du montage, et c’est là une raison 


de plus, raison majeure, de faire appel à d’autres. 
cristaux. | 
Gypse. — Le gypse présente un clivage principal 


et des clivages secondaires. ‘Il est facile de cliver. 
des lames minces dont les faces parallèles au clivage L 
principal seront courbées au rayon voulu. Une lame | 


offrir à la transmission un système réticulaire interne 
intéressant, par exemple (001) dont la constante 


réticulaire est de 2,8 À; (101) pour lequel d — 2,4; 
ou d’autres; ces plans sont normaux aux plans (010) 
des clivages principaux. Ceux-ci travaillent dans le. 
cas de la réflexion superficielle; la constante réti- 
culaire correspondante a été déterminée avec préci- 
sion, ainsi que sa variation apparenté d’un ordre 
de réflexion au suivant, par suite de la réfraction (°); 
elle est égale à 7,58 À. Les zones de longueurs d’onde 
qu’elle permet d'analyser avec le montage que je 
décris sont INCIQUÉeS dans le tableau. 


Quarlz. — Le quartz vi être taillé optiquement; 
il le peut être de bien des manières. H. Hulubei a 
indiqué, dans les lames ayant subi la taille de 
P. Curie, les caractères intéressants du système 
réticulaire interne (1340) : pouvoir réflecteur excel- 
lent en premier ordre, très faible en second ordre; 
constante réticulaire petite : d = 1,1776 À, incli- 
naison relativement faible : & = 13° 54. De nombreux 
autres plans sont utilisables dans la transmission, 


par exemple (1120), normal aux faces naturelles du 
prisme, avec. 2d—a—4,9g À: Avec une même 
lame ainsi taillée, les spectres peuvent se former par 


réflexion sur la face naturelle (0110) pour laquelle 
d — 4,24 À. : Cette constante a d’ailleurs été 
déterminée .avec une grande précision, ainsi que 
l'effet d'indice, ce qui est précieux (®). Les gammes 
de longueurs d’onde accessibles, dans notre montage 
(cas R), grâce à ce système réticulaire, sont spéci- 
fiées dans le tableau. 


Utilisation du spectrographe. 
graphe décrit pourrait être employé à l’analyse des 
rayons durs. Mais il n’est avantageux de travailler 4 
sous vide que lorsque l’absorption par l’air devient . 
importante. Ceci se produit déjà pour des lon- 


() On se reportera au traité de M. SIEGBAHN, Spektroskopie 
der Rôntgenstrahlen, 1931. 


VERSEL 5. 


Dispersions A) et AE ( Amin pour 100: 


“RAYONS X 


POUR LES 


Ras Dour B10AU En 


AE nÀ (4.X). - AA US) AE (eV). Méthode. 
EUX 6 MT TR EU 20,8 DO ï 
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gueurs one qui en de RO être ana- 
lysées par transmission, étant donnée la valeur des 
angles de Braggo correspondants. Aïnsi, je rappellerai 
que des rayons de 2 À sont absorbés de moitié 
= par 26cm d'air, ceux de 2,5 À par i4cm d'a 


cette zone de. longueur d'onde s’analyse bien par 
Dons à travers dés cristaux minces. À partir 
He A par: exemple, c'est-à-dire, déjà, pour les 


pectrographe. Fes 
Le tube est placé en ou T. Un calcul simple 


à étudier À = nÀ, réelle si 
n =: où «apparente » sin 1. Il faut donner au 
spectrographe la position d'ensemble © + 2, pour 
recueillir le spectre sur la plaque tone bheae 
placée à la position 2 Me 0) sur le cylindre de 
focalisation. 

Le, tableau donne ae exemples de disper- 
ons en longueurs d'onde (u.x.) et en énergies 
CV), pour quelques systèmes réticulaires commodes. 


M Lorsque la position du porte-plaque s’écarte 
: beaucoup de 90° (soit ox F)s les défauts 


‘de focalisation commencent à se faire sentir. Il 
: devient préférable de passer au cas R. Le tube est 


gueur d’ onde moyenne à 


faire que Sur. la surface concave du cristal; il n "Na 
_ donc qu’un système réticulaire utilisable par lame 
_ et c’est, pour le mica et le gypse, la face de clivage 
principal. Le spectrographe doit être incliné de 


ox suivant l'angle de Bragg 9n, calculé pour 


et ceux de 3 À par 8,5 cm. Or, une bonne partie de 


onne la valeur, de l'angle de Bragg, © pour la Jon- 


en u. x. et en électronvolts pour le mica et le gypse 


‘émissions K du cuivre, il ÿ aura nie à pomper le 


. tures convenables. L’élimination de ces rayonnements. 


là longueur d’onde moyenne. Le porte-plaque est 
placé à 2 (on) par rapport au centre de cour- | 


bure du cristal. Le tableau indique les déco 


jusqu'à la longueur d'onde limite, nas  COrres- 
pondant à l'angle limite de ce montage (81° 30). 


._ Notons qu’afin d'éliminer autant que possible les : 
rayonnements parasites, il y a lieu de disposer, 
entre le tube et le cristal, ainsi qu'entre le cristal 


et la plaque, des écrans métalliques percés d’ouver- 


est beaucoup plus difficile dans le cas, R que dans 


[cast 
Le pouvoir. séparateur est tel qu’ avec un bon mica 


les raies KP, et 6, du molybdène, par exemple, 


séparées par environ 0,5 u. X., 
une dispersion de 20 u. x. fmm. Il faut remarquer 


la dispersion dans l'échelle des énergies augmente 


placé en R. La réflexion ne peut, pratiquement, se 


que la dispersion des longueurs d'onde est beau- 
coup plus petite, en pratique, dans Ja région des 
grandes longueurs d’onde; mais que, malgré cela, 


beaucoup avec À Ce fait entraîne rapidement une 
grande. amélioration du pouvoir séparateur vers 
l’ultraviolet extrême, en ce qui concerne les sépa- 
rations énergétiques. Cependant, les émissions spec- 
trales prennent, aux faibles fréquences, l'aspect 
connu de bandes complexes dont la largeur intrin- 


sèque fait tomber le pouvoir séparateur apparent 


dans l’échelle des longueurs d’onde. 

Notre montage permet déjà, malgré le diamètre 
relativement faible du spectrographe, et donc dans 
de bonnes conditions de luminosité, l'analyse des 
niveaux de faible énergie, tels que les niveaux de 
conductibilité ou les niveaux réticulaires de Bril- 


sont dédoublées sous. 


= 


. d'absorption de. fréquences moyennes et dans la 


) région des rayons X mous. 


On notera, sur: les clichés de la planche I, la 
qualité du cristal analyseur (que traduit la finesse 
des spectres de rayons durs) et la largeur propre des 


- 1479 1985; 


185. 


_louin, qui se manifestent à la fois dans les spectres 
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au système HRvaes Aéeuons paramétriques… 


: H justi fie cette définition en montrant : 


AT NE de 


conséquences de leurs propriétés cinématiques. 


Introduction. — Depuis les travaux de Van der 
Pol, l'attention des physiciens s’est portée avec un 
intérêt croissant sur les phénomènes périodiques 


auxquels cet auteur a donné le nom d’ oscillations 


de relaxation. 


… Ces phénomènes périodiques sont caractérisés 
par une dissymétrie entre les phases de croissance . 


et de décroissance de la variable oscillatoire, la 
croissance étant lente et la décroissance rapide, ou 


réciproquement, et par le fait que le passage de la. 


phase de croissance à la phase de décroissance a lieu 
lorsque la variable oscillatoire atteint une valeur 
déterminée, qui constitue dès lors son maximum et 
définit l’amplitude de l’oscillation. 

C'est ainsi que varient. par exemple, le niveau 
de l’eau dans un vase de tantale alimenté par un 


mince filet liquide et la tension électriqüe d’un 


| DÉFINITION GINÉMATIQUE DES OSGILLATIONS Le RELA 
Par JEAN ABELÉ. 


Sommaire. — L'auteur propose une définition cinématique Fi oscillations de relaxation conduisant \ 


.3 = cosû — f(sin0), 


{désignant un: fonction arbitraire assu sais à Vérifier certaines conditions énoncées au cours de l'article. AE 


10 Qu’ elle généralise la ee de l’oscillation sinusoïdale comme projection, du mouvement cireu- D 


20 Que les to qui en résultent permettent, en Fire pour f une fonction: ARE un 
coefficient arbitraire, d’obtenir des oscillations évoluant, suivant la valeur de ce coefficient, de la Ai 
forme sinusoïdale aux formes les plus caractéristiques des oscillations de. relaxation ; ï 

30 Que ces équations conduisent à des équations dif'érentielles analogues à celle de Van der Pol, 
plusieurs propriétés dynamiques des oscillations de relaxation apparaissant dès lors comme es 


Vin, ME TETE 
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re 
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condensateur chargé à. travers une résistance. é Y 
et déchargé par une lampe au néon. LAN 

Une pareille forme d’oscillation suppose la variati 1 
à tout instant de l'énergie totale M l'oscillateur. 


 . ouvert sur l'extérieur, ce qui ne un 
dépendance de la vitesse de croissance ou de décrois- 
sance, et par suite de Ja période, vis-à-vis des gran- 
deurs caractéristiques du milieu où FOSCIaIQUE 
puise et restitue son énergie. NE MA Le 
Ainsi l'amplitude des oscillations de ibn est 
bien déterminée et constante, tandis que leur période 
est plus ou moins arbitraire et sujette à variations. 
On sait que l'inverse est vrai pour les oscillations 
sinusoïdales conservatrices d’énergie. . 
Malgré les caractères opposés de ces deux formes 
d’oscillations, il est possible d'insérer entre les types 


: 
La 


* 


+ 


Let 


f  discontinu de la phase de crois” 
ance à la phase de décroissance, toute une série 
e types intermédiaires qui les relient. Le grand 
iérite de Van der Pol:a été de montrer qu’une 
même Fe différentielle contenant un Ras 


HAE HIT + +h(D eo (1 


no sont des variables sans dimensions, où 
riables réduites: A gr 
st un coefficient positif, lui aussi sans dimension ; 


PE 


1(2) est une fonction paire, comprise entre zéro et 


croissant rapidement avec ce module lorsqu'elle 


D does à Vunité; ER 


nes € DAT V@TS ZÉTO. 
En posant pour simplifier a HE, 


" fé 


En n on obtient LAPeunt suivante 


De nee 


qui est la forme Mira et en classique de 
"équation de Van der Pol. Nous désignerons 


Van der Pol. 

; Bien avant les travaux de Van der Pol, Lord Ray- 
leigh avait étudié une équation différentielle non 
linéaire semblable, que nous écrirons d'une façon 


ne ARRET dz ]) dr ! px 
eee lun ]haten neo fe 


4 us de Lord Rayleigh. Une fonction de Me 


emplace la fonction de z présente dans le Hein 


_ conditions. 
Même sous sa forme simplifiée, la discussion de 
| Jéquation de Van der Pol exige des calculs d’approxi- 
mation assez longs ct des intégrations graphiqu:s 
pénibles dont les résultats, bien que très intéressants, 
sont limités. Si l’on en cherche la raison, on est tenté 
de répondre qu’elle se trouve dans le caractère 
non linéaire de l’équation, dont la solution est une 
fonction transcendante. 
. Cette réponse, en partie exacte, est en partie 
aire car la fonction sinusoïdale est, elle aussi, 
une transcendante. Mais celte transcendante nous 
_ est familière et nous disposons, pour la définir, de 
multiples ressources autres que l'équation différen- 


tre part l'oscil- 


ité pour les faibles valeurs du module de z, 


cinématique à la dynamique en formant, dans le : 


l'équation généralisée de Lord Rayleigh. 


‘équation (1) sous E nom d’équation généralisée de. 


par M définit une variable réduite égale au quotient 
lu second terme de (1) et est soumise aux mêmes . 


tielle du ie DRE à siens constants. Elle 
s'introduit en particulier en mécanique , comme 
_ projection du mouvement circulaire uniforme, c’est- 
à-dire par une définition purement cinématique. É 
L'absence d’une semblable définition des oscil- 
lations de relaxation nous paraît être la principale 
lacune de nos connaissances en leur domaine. Le 
but du présent travail est de la combler. 
Nous montrerons, dans une première partie, que 27 
la projection sur un axe fixe parallèlement à un ate 
mobile d’un mouvement circulaire non ArOrRIS 
convenablement, choisi, permet d’engendrer des 
oscillations possédant les propriétés ci-dessus décrites 
comme caractéristiques des oscillations de relaxation, 
et nous déduirons de ce mode de projection des 
formules générales qui en de une expression 
analytique. 
Au cours de la seconde partie, nous construirons, 
dans un cas particulier choisi comme exemple, les 
courbes correspondant aux différentes valeurs d’un 
paramètre variable, et nous donnerons la formule 
permettant de calculer les fréquences ainsi que la 
courbe correspondant à cette formule. gi 
Nous passerons, dans une troisième partie, de He 


même cas particulier, l’équation différentielle du 
second ordre, et nous vérifierons qu’elle se ramène, 
soit à l'équation généralisée de Van der Pol, soit à 


\ 
Î 


I. Génération cinématique des oscillations de 
relaxation par projection sur un axe fixe, 
parallèlement à un axe mobile, d'un mouve- 
ment circulaire non uniforme. 


1. Conventions Lun au mode de projec- 
tion. — (Considérons, comme pour la définition 
des fonctions trigonométriques, un vecteur uni- 
taire OM décrivant l’angle 0 à partir de OZ dans 
le sens indiqué par la flèche (fig. 1). Sa projec- 
tion Om — z sur le diamètre ZZ' du cercle engendré 


par une amplitude À de la variable non réduite qui 
lui correspond. Cette dernière sera. désignée par le 
produit (A2). 

L’axe YY', parallèlement auquel est effectuée 
la projection, peut osciller autour de O. L'angle 8 qu'il 
forme avec ZZ' est donc variable, ainsi que son 
complément 9. 

Les coordonnées z et y sont données par les 
formules 


sin À 


Cr ALL AP 
| çosa 


RU ER, pee en 


sin f5 


coS(0 + a) 
WCOoSa 


Sin 
- sinf 


(3). 


à 


Il suffit d’adjoindre à l'expression donnant z 
une relation entre un des angles x ou 5 et Ô pour 
que z soit complètement déterminé en fonction de 0 CR 


x 
we 


98. ” 


L'oscillation de l’axe YY' peut être apériodique 
ou périodique et, si elle est périodique, elle peut 
être asynchrone ou synchrone de la rotation du 
vecteur OM; nous avons utilisé ailleurs (?) la première 
de ces trois hypothèses, nous: n’envisagerons ici 
que la dernière. ! 


Ai 


« RER 
Fig, 1. — La lettre $ désigne l’angle ZOY, et non MOY. 


2. Génératrice d’une oscillation de relaxation. 


— Lorsque les deux mouvements sont synchrones,. 


il est avantageux de déterminer z en fonction de 0 
par une relation entre le vecteur Op = y et un des 
angles « ou B. Cette relation est l'équation en 
coordonnées polaires de la courbe C décrite par le 
point p, projection de M sur l’axe oscillant. 

Si de l'équation de C en coordonnées polaires nous 
passons à son équation en coordonnées cartésiennes 
rectangulaires dans le système d’axes (OX, O7) 
et si nous désignons par f(x) l’ordonnée qp, la figure r 
montre que la courbe C jouit d’une propriété remar- 
quable. On a en effet 


:=pM=qM—f(z), 
gM = cos8, æ = sind, 
3 = cos — f(sinb). (4) 


Construisons dans un système d’axes rectan- 


gulaires (0, 0, z) la courbe :d’équation (4) et la 


sinusoïde z — cos 0 que nous appellerons sinusoïde 
de référence (fig. 2). La fonction (sin 0) est égale 
à la différence des ordonnées de ces deux courbes 


() Thèse, Paris, 1943 et Annales françaises de Chrono- 
métrie, 1942, fasc. ‘3-4; 1re Partie, Chapitre V, L’objet de 
cette thèse est de démontrer la possibilité de réaliser, en théorie 
et en pratique, un oscillateur à la fois non linéaire et rigou- 
reusement sinusoïdal, susceptible en conséquence de concilier 
la stabilité de l’amplitude et la stabilité de la fréquence. 
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auquel à lieu Ja projection. 


trice en vue d'obtenir des oscillations douées des. 


mènes de relaxation. : 


ape" /! RE SAC 


' 


À 3 y te A Ar 4! h 
(indiquée avéc son sens par des flèches). Elle mesure … 


la transformation subie par le mouvement projeté 
de M, du fait de l’oscillation de l’axe parallèlement 

À la courbe C, définie par la fonction f(x) qui repré-. 
sente son ordonnée dans le système d’axes (OX, OZ), 
nous donnons le nom de génératrice de l’oscillation 
de relaxation. ‘ te 


3. Conditions imposées à la génératrice. — 


Déterminons les conditions à imposer à la généra- 


propriétés essentielles des oscillations de relaxation. 


La première de ces conditions est que la fonc- … 


tion f(x) et ses deux premières dérivées s’annulent 
en même temps que x, en sorte que la génératrice 
passe par l’origine et y ait un point d’inflexion à 
tangente horizontale. En vertu de cette condition, 
les maxima et minima de la courbe d’équation (4) 
sont égaux à ceux de la sinusoïde de référence et 
ont lieu pour les mêmes valeurs de 0. De plus, les. 
deux courbes se traversent à chaque maximum 


et minimum (cf. fig. 2), en sorte que, si l’oscillation : 


de relaxation éroît plus lentement que la sinusoïde  * 


immédiatement avant un maximum, elle décroîtt … 


plus rapidement immédiatement après, et réci- … 


proquement. 


Ainsi sont sauvegardées l’invariance de l’ampli- 


tude et la dissymétrie de croissance, propriétés que 
nous regardons comme caractéristiques des phéno- 


La seconde condition concerne l’évolution de la 
génératrice de part et d'autre de l’origine et en 
restreint le domaine à l’espace limité par deux demi-. 
cercles tracés de Z et Z' comme centres avec l’unité 
comme rayon et tangents en O (fig. 3). 

En vertu de cette condition le module de z ne 


peut être supérieur à l'unité et les maxima et minima 


relatifs atteints pour 0 — kr sont aussi des maxima 
et minima absolus. 

De plus, la partie de la génératrice effectivement 
parcourue par la projection de M sur l’axe mobile 
est limitée par deux droites parallèles à ZZ' passant 
par X et X’, en sorte que les points où la génératrice 


MEET ACTE. SON 


# 


4 


À ; : te. 
Le) 49 $ Fe à 


oupe ces droites ‘sont, pour le point mobile p. 


des points de rebroussement. 

Nous imposerons enfin, mais uniquement en vue 
de simplifier notre étude, une froisième condition 
- à la génératrice, en exigeant qu’elle admette l’origine 
comme centre de symétrie. Il y a alors dans le 
_ plan (z, 0) symétrie par rapport à l'axe des 0 des 


deux moitiés, positive et négative, de la courbe 


“ représentative de l'oscillation si on les déphase 
_ d’une demi-période (cf. fig. 2). 


EN 


NN 


< Fig. 3. 


Une conséquence de cette symétrie est que la 
période de l’oscillation de l’axe est égale à la moitié 


_ de la période de rotation du vecteur OM. 


4. Loi de mouvement du vecteur tournant. — 


Soit & la vitesse angulaire constante du vecteur OM 
_ lorsque, « restant constamment nul, la projection 


de M sur ZZ' se meut d’un mouvement sinusoïdal. 
Nous adopterons le produit (sf) comme variable 


_ temporelle réduite. On a alors 


dô A 
dot) 


© 


et (ot)= 0 — 0. ( 


Recherchons la loi de mouvement qu’il convient 
d'adopter lorsque l’axe YY' est animé d’un mou- 
vement oscillatoire synchrone du mouvement de 
rotation. 

Une question préalable se pose. Nous avons vu 
que l’oscillation de l’axe YY” transformait, moyen- 
nant quelques conditions simples imposées à la 
génératrice, la sinusoïde de référence en courbes 
présentant les caractères des oscillations de relaxa- 
tion. Puisque le mouvement uniforme défini par 
les équations (5) conserve la forme de ces courbes 


_ lorsqu'on passe de la variable 0 à la variable (w1), 
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il peut sembler inutile d'introduire entre ces variables 
une relation non linéaire. Deux raisons cependant 
concourent à rendre cette introduction nécessaire. 

La première raison est tirée du caractère de varia- 
bilité de la période que nous avons signalé dans notre 
description initiale des oscillations de relaxation. Les 


facteurs concrets responsables de cette variabilité 


sont représentés mathématiquement par les para- 
mètres qui figurent dans la fonction f(x). La loi 
de mouvement du vecteur tournant doit donc 
dépendre de cette fonction. 

La seconde raison est tirée de l’impuissance du 
système de projection par axe oscillant à rendre 
compte, à lui seul, des formes extrêmes que peuvent 
revêtir les oscillations de relaxation. 

Pour démontrer cette impuissance, donnons à la 
génératrice la limite des formes qu’elle peut recevoir 
en la constituant par les deux quarts de cercle ACO 
et ODB et étudions les variations de l’élongation 
lorsque l’extrémité M du vecteur tournant décrit 
à vitesse constante l’arc de cercle X’ZD (fig. 3). 


L'élongation étant égale au segment PM mené 


parallèlement à OZ de la génératrice au cercle, on 
voit qu’elle reste constamment égale à l’unité pendant 
que M décrit l’arc X'Z et qu’elle décroît de l’unité 
à zéro pendant que M décrit l’arc ZD. Tandis que 


la croissance infiniment ralentie qui précède le. 


passage de M par Z donne bien la limite du ralen- 
tissement que l’on observe sur les parties à faible 
pente des courbes obtenues par Van der Pol, la 


. vitesse de décroissance qui suit le maximum, loin 


d’être une limite supérieure, est inférieure à celle 
des parties à forte pente de ces mêmes courbes et 


l’on ne comprend pas qu’une discontinuité puisse 
‘séparer ces deux phases de l’oscillation. 


Ces caractères ne peuvent être obtenus que par 
une variation de la vitesse angulaire du vecteur 
tournant. Cette vitesse doit passer, plus ou moins 
rapidement, d’une valeur inférieure à l'unité à une 
valeur supérieure lorsque M passe par Z dans le 
sens indiqué sur les figures 1 et 3. < 

L'hypothèse la plus simple permettant d’obtenir 
ce résultat consiste à donner à la vitesse angu- 
LT. la valeur que possède le rapport de la 
pente de l'oscillation de relaxation à la pente de la 
sinusoide de référence dans le système de réfé- 
rence (0, z, 0) de la figure 2. 

Nous poserons donc | 


laire 


dz 
dë d8 
dtor) CET sin (6) 
d’où 
_ sin8 .. 
dot) = _. di. (Gbis) 
d5 


Cette convention recevra sa pleine justification 
de l'étude des Cas particuliers. 


5. Équitiue pe CR UraIes des 
! oscillations de relaxation. -— L’adjonction à 
_  léquation (4) de l'intégrale de l’équation (6 bis) 
constitue un système d'équations paramétriques 
reliant l’élongation z à la variable temporelle 
réduite ot 


z = cosô — (sin), A (4) 
o Ê 
— sinÔ 5 
= f Te à. (7) 


Ces équations résultent de la définition cinéma- 
tique que nous avons donnée des oscillations de 

relaxation. Mais on pourrait les poser «a priori, 
ip dans l'esprit de la méthode axiomatique, avec les 
_ conditions précédemment énoncées pour la fonction 
de arbitraire qui y figure, conditions auxquelles il 
faut ajouter celle d’après laquelle l'intégrale du 


toujours être positive et finie. On aurait ainsi une 
définition analytique formelle des oscillations de 
relaxation: les considérations qui précèdent n’en 
donneraient qu'une représentation. Au lecteur de 
‘ choisir le point de vue qu’il préfère (?). 

Nous allons appliquer les équations (4)- (7) ya 
un cas particulier, en étudiant successivement 
_ quelques propriétés cinématiques et dynamiques de 
- Ja famille d’oscillations obtenue en adoptant pour 
de la génératrice l'équation la plus simple satistaisant 
aux conditions imposées. 

| ; \ 


II. Propriétés cinématiques 
_de la famille d’oscillations définie 
| par la is ot A d'équation f(x) — a. 


D Us Équations paramétriques. —— La fonction f(x) 
la plus simple qui s’annule ainsi que ses deux pre- 
_mières dérivées en même temps que x et qui soit 
_ symétrique par rapport à l’origine est ax. Les autres 
j conditions pourront être satisfaites en limitant le 
domaine des valeurs du module de a. 

Les équations (4) et (7) deviennent alors 


3 = cos — & sin*0, (ei 
0 FRA) 
d0 
ed eus x ) 
(wt) dé 1 a sin Ü cosô À 


() I se trouvera peut-être des lecteurs qui estimeront 
arbitraire le choix des propriétés (constance de l’amplitude, 
variabilité de la période, dissymétrie de croissance et de 
décroissance) par lesquelles nous avons caractérisé les oscil- 
lations de relaxation; choix qui est au point de départ de 
notre construction cinématique. 
dans les Principes fondamentaux de physique théorique de 
J. L. Desroucuess, t. I, Chap. V, les observations de cet auteur 
sur la nécessité préalable à toute axiomatique d’une Synthèse 
inductive (cf. en particulier les paragraphes 10 et 11). Nous 
avouerons ensuite sans aucune gêne que nous ne prétendons 
point convaincre le lecteur, mais que nous espérons entraîner 
son adhésion en lui permettant de saisir intuitivement l’idée 
générale qui nous à guidé, 


second membre de (7) prise entre o et 27 doit 


à a —o. Si nous portons en ordonnée, avec une 


Nous les invitons à lire 


=— arctg 


+] 
En faisant 0 — Set 6, =— <et en doublant le. 
résultat du calcul on obtient la période GT). pe 


OT) = + 


Désignons par ÆF, la fréquence correspondant 


à ; F 
échelle graduée de o à x, les fréquences relalives —- Pr 
0 


et si nous portons en abscisses | a| avec une échelle: 
dilatée par rapport à la précédente, le coefficient 
de multiplication des intervalles étant 3/2, la courbe 
donnant la variation des fréquences cst un quart de 
cercle de centre O (fig. 4). 


3. Variation de forme des courbes en fonction 
du coefficient a. — Pour comparer les courbes 
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III. Propriétés dynamiques 


t a, il est intéres > changer e de _ de la famille d’oscillations définie 
1bsci: e façon à r es périodes identiques. par la génératrice d'équation f(x) — axi. 


Il suffit pour cela de substituer à la variable (ot) 


e rapport 7 dont l'expression s'obtient en divisant  Î- Expressions générales de la vitesse et de 
CR DE PRG D ne ur l'accélération. — De l'équation (6) qui détermine | 
bre à membre les équations (10) et (11) an ; À nu 

je NN la vitesse angulaire —— du vecteur tournant OM Fa 


d(wt) 


on déduit | PE 
CAN AT RS 


2. Formation de l'équation difiérentielle du 

_ deuxième ordre admettant (8)-(9) comme solu- 

tion périodique. — Nous allons appliquer les 

. formules générales (14) et (15) à l'équation (8) 

T = cos 0 — a sin*Ü et nous montrerons que la rela: 

Ê TT AS d?z | 45 

Mo ta CO ir 

le signe de a, sous la forme, soit de l'équation gén 

ralisée de Van der Pol, soit de l'équation généralis 

de Lord Rayleigh RE A ue 
dz ie d 

d(wt) dur) 

= sinÜ(1 + 3asin0 cos0 }?, 


Y 


peut être mise, suivan 


(cosi — & sins 6). 


 d?z 
d(orR à re 
nl + 6a sin0(1+ 3asin8 cos 0)(cos?0 — sin?0)]. 
X (14 3a sind cos®) | FUX 

=— [cosÙ +9asin cos 0 — Gasin*0] 
X (14 3asinb cos0 }? PRE SENS 
= — 3a(2sin0— 3 cos®0)[—sin0(1+3asin8 cosô)?] 

— cos0(r1 + 3 a sin0 cos}. PS IQS 


— | cosä(1 +3asin0 cos 0 )? 


1 


… L'’équation différentielle cherchée est donc, 0 pouvant 
être considéré comme une fonction, soit de z, soit . 


dz | ne 
PTIT ous pie aa) , 1 cp 
d?z | : dz: 
——— É 9\S1N20 2" 2 
dur} + 3a(2sin 3 cos Mara M 
+ cosû(i + 3a sin 0 cos0 2 = 0. (17) 


“TRS ‘a Supposons d’abord a <o, nous écrirons (17) . 
Nous poseronsf0, — o, ce qui revient à prendre sous la forme suivante : ; k * 
pour origine l’abscisse du maximum de z. Il sufit 4. RENTE : à 

alors de faire varier 6 de Ô à x pour obtenir la rs +balie, cost) 

moitié de la courbe, l'équation (13) donnant les | 400801 38m 0 cosb y 2 0. (AT UE 
_ abscisses et l’équation (8) les ordonnées. L'autre RTE 


moitié s'obtient par glissement et retournement. Il suffit alors de poser 
Les trois courbes de la figure 5 ont été obtenues Hd & U18) 
en donnant à a les valeurs + 0,2, + 0,4, + 0,6. di 
En adoptant les valeurs négatives correspondantes dés) 2 6 CE, (49) 
on obtiendrait des courbes symétriques des précé- ; 2 ou 
dentes .par rapport à Oz. cos0(1 + 3 @ sin0 cos )? ji | 
! L’analogie de ces courbes avec celles qu’a obtenues De (: a : sint cost) Mean (20) 


Van der Pol est frappante. 


MNT Ta NP Ile 
pts » | ot 


402 0 Port Us 


pour que l'équation (17 bis) se. ramène à l'équation 


généralisée de Van der Pol : 


d?3 
d(wt)? 


+3) = 0. 
on PHARE à (0) 


= (2) 


Supposons maintenant a = o, nous écrirons (17) 
sous la forme suivante : 


Res Ô in? 6 2 
d(wiÿ —oefi— ge d'wt) 

+ cos8(1+— 3asin0cos0}?=o. (174) 
Il suffit alors de poser 
(e> 0), (21) 
5 dz 
3 SN 0 = ge fe 29 
ON NE ni ri 


cosÜ(1 + 3asin0 cos0 )? 


94 — € 


= cos (+ ;sino cost)" = Vo(z) 


T 


(33) 


pour que l'équation (17 fer) se ramène à l'équation 
généralisée de Lord Rayleigh : 


one els HAE dz dz DAS 270 
PAR e [élan + Hits) 20 (2) 


3. Variation des fonctions v,(z) et d.(2). — En 
associant l'équation (8) aux équations (19) et (20) 


Ait. À 
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nous obtenons un système d'équations paramé- 
triques qui permet de construire les courbes repré- 


sentant les variations des fonctions o,(7) et (2). 


Z —= Cos0 — a sin30, (8) 


1 
p1(z) = = cost, 


(19) 
Y1(2) = cos0(1 + 3 a sine cos 0 }?. (205 
La même méthode permettrait d'étudier les fonc- 
: dz 
tions ®» | PTE res) et d (2) à l’aide des équations (8), (16), 
(22) et (23). 
La figure 6 donne deux courbes représentant la 


PRET A (1° AN OC Er EL FRE Put RS Te INPACT TA ., PEN 1) ES 
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variation de 10 : 
rompu), 22 Pour 4 ——0,2 (en trait continu). LA 


première est une parabole à axe vertical dirigé vers. k 
le haut, dont le sommet est à l'origine des coor- 
données. Elle est identique, à une constante près, ! 
à la parabole représentant la fonction À? qui figure 
dans le coefficient du second terme de l'équation 
simplifiée de Van der Pol. La seconde courbe résulte 
d’une déformation de la première, les extrémités A 


et C restant fixes et le sommet subissant le long de 
l'axe des z une translation égale. .à — a. 


LP 


Cette déformation ne modifie pas essentiellement 
le mécanisme de l’amorçage des oscillations et de la 


limitation de leur amplitude. Pour les faibles valeurs 


de |z|, w,(2) reste inférieur à l’unité et l’amortis- 


sement est négatif, Ce qui provoque l’amorçage 


des oscillations. Pour les valeurs élevées de |z|, 
o.(z) est supérieur à l’unité et croît rapidement, 
l’amortissement est alors positif et rapidement 


croissant; de là résultent la limitation de l’amplitude 


et l'existence d’un régime permanent. 

Remarquons que, du fait même de notre définition 
des oscillations de relaxation, l’amplitude du régime 
permanent est rigoureusement constante, alors 
qu'elle est légèrement variable dans l’oscillateur 
régi par l'équation simplifiée de Van der Pol lorsque 
l’on fait varier &. 

La figure 7 donne de même deux courbes repré- 
sentant la variation de d,(z) : 1° pour a — 0 (en 
trait interrompu), 2° pour a — —0,2 (en trait continu). 
La première est une droite passant par l’origine 


10 na a=0 1 nr nters g 


simplifiée de Van der Pol. La seconde peut être 
considérée comme une déformation de la première, 
la partie médiane subissant, avec une légère ondu- 
lation, une translation vers le bas et les parties 
_ extrêmes étant fortement incurvées. Cette défor- 
: mation appelle deux remarques. 


1° On sait que dans le cas de l'équation linéaire 
_ à coefficients constants le facteur d'amortissement, 
s'il est faible, 


terme. De même, dans le cas de l'équation non 
linéaire, le troisième terme joue un rôle plus impor- 
tant que le second dans la détermination de la 
période. Dans l'équation simplifiée de Van der Pol, 
qui respecte la linéarité du troisième terme, la 
à période varie peu pour les faibles valeurs de € et 
elle ne tend vers l'infini que lorsque e tend lui-même 
_ vers l’infini Dans l'équation généralisée (17 bis) 
Ja non-linéarité du troisième terme entraîne une 


fonction de € et celle-ci tend vers l'infini quand € 
tend vers 4. : 


20 L’altération progressive du troisième terme 
pour des valeurs croissantes de € suit normalement 
celle du second terme dans l’équation de tout oscil- 
lateur concret entretenu par l'introduction d’une 
résistance négative. Une semblable résistance n’est 
en eflet presque jamais une résistance pure, c'est- 
-à-dire que la force d'entretien n’es{ presque jamais 
rigoureusement en phase avec la vitesse. 

Elle se comporte comme la somme de deux compo- 
_ santes, l’une en phase et l’autre en quadrature avec 
la vitesse. La première rend non linéaire le second 
terme de l’équation et la seconde rend non linéaire 
le troisième terme. | 

La composante en quadrature peut être diminuée 
ou partiellement compensée par, divers” procédés, 
elle n’est jamais complètement négligeable; elle 
l'est d'autant moins que la valeur moyenne de la 
force d'entretien est plus élevée, c’est-à-dire que € 
va-en Croissant. Ainsi l’équation que nous avons 
proposée dans l’introduction sous le nom d’équation 
généralisée de Van der Pol et qui se distingue de 


1 


ion 4 ee ie dans Le Lrisiste En de l ie 


influe peu sur la période qui est. 
surtout déterminée par le coefficient du troisième 


_ variation beaucoup plus rapide de la période en. 


Péquation simplifiée du même auteur par la Cabot. 
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#4 


_tution de fonction d,(z) à z, est-elle plus proche 


des conditions réelles d’entretien. Il est d’autant 
plus intéressant de constater que c’est à une semblable 
équation que conduit, dans le cas particulier étudié 
ci-dessus, la définition des oscillations de relaxation 


que nous proposons. 


Conclusions. — Bien que nous n’ayons étudié 
dans la [feet la IIIe Partie de ce travail qu’un cas 
particulier et que notre intention soit de généraliser 
cette étude avant de formuler des conclusions 
d'ensemble, nous croyons cependant que les résultats 
déjà obtenus jettent un jour nouveau sur la méthode 
avec laquelle ïil convient d’aborder l'étude des 
phénomènes de relaxation. 

S'il s’agit de phénomènes concrets, tels que les 
battements du cœur, l’évolution statistique de 
certaines espèces biologiques, etc., les équations 
auxquelles nous avons abouti permettront de 
rechercher la forme de génératrice qui leur corres- 
pond, et l'exemple étudié ci-dessus nous autorise à 
penser que des génératrices relativement simples 
aideront à représenter, à comparer et à classer ces 
phénomènes. 

S'il s’agit d'étudier les propriétés théoriques des 
oscillateurs entretenus, la création d’une cinématique 
des oscillations de relaxation permettra de classer 
ces propriétés en deux catégories : celles qui pro- 
viennent de la non-sinusoidalité des intégrales et 
qu'il nous apparaît désormais possible de déduire 
directement du système d’équations que nous 
proposons, et celles qui proviennent de la non-linéarité 
des équations différentielles et qui exigeront toujours 
une étude de ces équations (c’est le cas par exemple 
de la synchronisation et de la démultiplication de 


fréquence) (). Mais même pour l'étude de ces 


dernières, une cinématique des oscillations de 
relaxation sera utile et permettra en particulier de 
mieux comprendre le mode de formation des équa- 
tions différentielles. 


(5) Le fondement de la distinction de ces deux catégories 
de propriétés repose sur l’existence d'équations non linéaires 
admettant une solution de régime permanent parfaitement 
sinusoïdale, les oscillateurs régis par ces équations étant 
susceptibles d’être synchronisés et démultipliés. Cf. Y. RocArb, 
Dynamique générale des vibrations, Ghap. XV, Paris, Masson, 
1943 et notre thèse déjà citée. 


Manuscrit reçu le 24 mars 1945. 
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Ü NE ca conduit à une valeur de y = =. 


température réduite 


Th 


rA Potentiels d'interaction. — - La théorie ciné- 
FR tique des gaz représenta d’abord les molécules par 
_ des sphères rigides et élastiques de diamètre 9; 


d'interaction. E. de deux telles Mn lÉbeS dont les 


() 


de la théorie fut de prévoir l'indépendance du 
coefficient de viscosité d’un gaz avec la pression. 


racine carrée de la température absolue WT, alors 
que l’expérience indiquait une variation nettement 
plus rapide. 
Maxwell considéra alors les molécules comme 
_ des centres de force répulsive proportionnelle à une 
puissance négative de la pue (potentiel d’inter- 
action r b en r-*) et prit $ — 4, valeur pour laquelle 
les calculs se simplifient énormément; il trouva 
le coefficient de viscosité proportionnel à T, variation 
cette fois plus rapide que celle donnée par l’expé- 
rience. Chapman et Enskog entreprirent les calculs 


INTERACTIONS MOLÉCULAIRES ET THÉORIES DE LA VISCOSITÉ he 
ET DE LA DIFFUSION THERMIQUE DES GAZ TA Na AN Va 


Par ANDRÉ FOURNIER. ÿ 


Sommaire. — La meilleure représentation de l'énergie potentielle d'interaction moléculaire est la 
forme de Lennard-Jones. En adoptant les exposants 12 et 6, on peut faire l’étude des propriétés d’un 
gaz à l’aide de deux constantes, la distance r, à laquelle le potentiel s’annule etla profondeur E,, du trou 
de potentiel. Celle-ci est reliée à la température cr itique; les données expérimentales ne sont pas suffisantes 
pour la rattacher à l'énergie de vaporisation. Le simple potentiel répulsif en r-* prévoit un coefficient de ONE. 
: viscosité n proportionnel à une puissance Y y de la température absolue et un facteur de diffusion ther- FPS HE 
We ‘mique indépendant de la’ opese Le potentiel de Lennard-Jones, en accord Avec les résultats 


d In À 
. un facteur de diffusion ther A ble avec la température. Étant donné le nombre de constantes EE 
caractéristiques du gaz, on peut mettre les résultats sous une forme réduite où n'intervient que la UT A 


+ Gette forme de la Loi de états correspondants n est qu ’approchée ; celatientàa lé Fa & 


. ce que l'approximation de la forme ‘exponentielle de la partie répulsive du potentiel (valeur finie 
pour r — 0, troisième constante caractéristique) par une forme en r-* bras DrORANIENENE être faite 
avec un exposant s variant un peu d’un corps à l’autre. : : 


 l’énergie potentielle (ou, pour abréger, le potentiel). 


centres sont distants de r est nulle pour r > a et. 
‘a in finie pour r a (y (Ag. 1 a). Un des grands succès 


. Mais ce coefficient était trouvé proportionnel à la 


décroissant quand la température augmente et à 


4 2 


compliqués valables pour une valeur quelconque 
| es S examinèrent Fe di ue de XIE OA 


no L ù UE ; 

Dans tous ces éalculs: les nn Hensleeule es 
‘attractives étaient négligées; ces forces de cohésior 
sont Loeue du. terme de por interne de 


tantes dans les états ne ae ou que) 
Elles sont en première approximation négligeabl s 
aux distances intermoléculaires moyennes d'un g; 

Hans Les conditions ordinaires ; mais Sutherland 


pas de sans on etil en aise pour le cbeicien 
de viscosité une formule bien vérifiée dans d'assez 
grands ‘intervalles de température. Le potentiel | ke 
utilisé, avec partie attractive en 77", est de la 
forme 1 c) On peut d’ailleurs distinguer deux cas ere 
le potentiel de Sutherland pour lequel la partie 
attractive est très faible; la profondeur du trou 
de potentiel EŸ intervient seule dans les résultats 
dés calculs et non l’exposant s'; et le potentiel 
plus général de Langevin-Keesom pour lequel cet. 
‘exposant subsiste dans les résultats [2]. 

Quand le potentiel attractif devient de plus en 
plus important, on est conduit à prendre pour le. 
modèle 1 b un exposant s de plus en plus petit 
(age, à 5,6 pour Host et ce modèle repré- 


} 


; d et Lennard-Jones l'a représenté 
rmule du type : 

RE Nr. He. (1) 
interactions sont d’origine électrique, A 
s quand les couches électroniques des deux molé- 
arrivent . à ne l'une de l'autre, ét 


as de l Ééhcn: le terme répulsif est exponentiel, 
s la. forme en r* constitue une approximation 


voit que ce or à peut être considéré comme 
superposition de répartitions dipolaire, quadri- 
p. rèz _octopolaire, etc, d’ importance décroissante. 


entes (par exemple dipôles permanents des 
ales polaires), mais il sy ajoute en tout cas 


" de en UE plusieurs sortes de es attrac- 
ves [A]: = 


E npemients  — permanents = Se À 
S. nas d'interaction sont au premier ordre 
r-9+1#1) et au second ordre en r7#7*7+1), Entre 
ôles permanents (q — q' = 1), l'effet du premier 
re est nul; il subsiste l'effet du second ordre 
n r-6 dont lexpression est, en désignant par p 
es moments dipolaires et par k la constante de 
Boltzmann, F 


PL. (2) 


>our les molécules non polaires, le premier terme 
non nul en moyenne est le terme du second ordre 


rares, on peut affirmer que les moments quadri- 
_ polaires permanents sont nuls). Quand la tempé- 
_rature augmente, ces forces d'orientation diminuent, 
car l’agitation thermique gêne de plus en plus l'arrêt 
._ dans les positions d'attraction (alignement des 


“LRO 


j 20 Les forces d’induction de Debye entre arran- 
Gr gements polaires permanents et arrangements polaires 
- nu par les permanents. L”° Dee du potentiel 


art d molécules, A ss 


quadripoles permanents en 7-1 (pour les gaz: 
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7e attraction entre dipôle permanent et dipôle induit 
est Æ 
a p'? + à p? 
— 76 £ (@) 
où les & sont les polarisabilités. Ces forces d’induction 
sont indépendantes de l'agitation thermique, la 
répartition induite variant comme la répartition 
inductrice. 


30 Les forces de dispersion de London. Ce sont 
les seules pour les molécules ne portant pas de 
distribution polaire permanente, mais même pour 
les molécules polaires ce sont les forces prépon- 
dérantes. Leur origine est analogue à celle des forces 
précédentes, mais les dipôles inducteurs sont instan- 
tanés. La mécanique classique, ne prévoyant pas 
d'oscillation résiduelle du système des charges 
autour de la distribution d’équilibre non polaire, 
ne pouvait en rendre compte. Ces forces résultent 
en somme de la tendance qu'ont les charges de 2 molé- 
cules à effectuer leurs mouvements en phase. Le 
potentiel d'attraction est donné avec une bonne 
approximation par l’expression 


3 ax. V;V; : 
e ee (4) 
D TS V; + Æ 
où les V; sont les potentiels d’ionisation, suivie 
comme précédemment d’une série de termes 


enr T0 


En résumé, dans tous les cas le terme principal. 
du potentiel attractif est en r—$. 


3. Équation d'état. —— Prenons le cas d’un gaz 
pur. Foutes les molécules sont identiques ou diffèrent 
tout au plus par les différents isotopes des atomes 
constituants. Les structures électroniques de deux 
composés isotopiques sont identiques, puisqu'elles 
dépendent essentiellement des charges nucléaires et 
non des petites variations de masse (sauf peut-être 
pour les tout premiers éléments). Les potentiels 
d'interaction entre molécules isotopiques sont donc. 
identiques. 


L’interaction moléculaire est représentée par un 
potentiel de Lennard-Jones avec $ = 12 et s' — 6. 
Exprimons-le sous une forme qui mette en évidence 
deux constantes caractéristiques : la valeur absolue 
du minimum d'énergie potentielle Æ,, profondeur 
du trou de potentiel, et la distance r, pour laquelle 
le potentiel s’annule, distance assimilable au diamètre 
de la molécule, étant donné la rapide montée du 
potentiel répulsif 


Ar} (En | (2) -() | 
 F(r)= ee p(e) (2) 


(à) 
La force 


(6) 
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s’annule pour r —7r, tel que 
RE AE A) 


Quels renseignements peut-on tirer de l'équation 
d'état sur la valeur de E,» ? Cherchons l'équation 
d'état valable pour de faibles écarts à la loi des 
gaz parfaits, donc pour de faibles pressions, c’est- 
à-dire l'expression de P V.en fonction des puissances 


: I : 5 F au 
croissantes de v° dont les coefficients s'appellent 


coefficients du viriel. À 
D’après le théorème du viriel de Clausius, l'énergie 
cinétique d’un gaz est égale à son viriel, valeur 


moyenne de — = Z(2X + yY + zZ), où la somme È 


est étendue à toutes les molécules du gaz, et où X, 


Y, Z sont les composantes de la force agissant sur 
chaque molécule. On en tire 


PV=NKT+ RS (r) (8) 


| (la somme double est étendue à toutes les paires 


de molécules; N est le nombre de molécules, nombre 
d'Avogadro si l’on prend pour V le volume molaire). 

En tenant compte de la probabilité de trouver 
#2 molécules à la distance r, l'équation précédente 


© devient 
A E(r) 


27 N ue es 

er Fr) e k1 | 
27 N. R DA $ 

ÆENE r} ne a [ 1? [, —e NT | | (9) 


On appelle point de Mariotte-Boyle Ty la tempé- 
rature pour laquelle le second coefficient du viriel 
s’annule. Dans le réseau d’isothermes du dia- 


PEN ur FE 


gramme P, PV, c’est la température à laquelle le 


lieu des minima de PV coupe cet axe. Lennard- 
Jones [5] a indiqué pour ce coefficient un dévelop- 


pement en série valable pour toute valeur de s 


et. de s”. Pour s = 12 et s:— 6, la condition de son 
Ba 
annulation est, en posant ET = %, 
PNR Zs iN æ 3\ æ2 5\ &° 
EE ES € A en LE y AREA A à Je + 
. (5) : (7) 2 G 3 
DM EI O2 UT ‘ 
+r( ; + (10) 
E 
soit + — 0,5/ ou, en posant 0, — T? 
(u 
0» En 
Tu A kTm tue sie) 


Quelle relation y a-t-il entre le point de Mariotte- 
Boyle et la température critique ? D’après l’équation 
réduite de Van der Waals, 


7 MINE 
LES 9 
12 Fe RAM ULe (ee) 


On en tire 
Rae 
dP" P'=0 Fe f | 12 
7 Fe > 
Ta —=V6 —= 2,40 


attractif 4E»r dont les constantes sont fournies 


. du viriel et covolume) avec les expressions du para- 


Pour les bee pressions, la meilleure 
réduite est celle dep: Berthelot- Clausius 


TOME ra 39 D 
re Den es re 
et ee à ae 


Les lettres primées désignent des quantités réduites. 


Les quelques valeurs a de varient entre 2 1 
et 2,8. Adoptons 2,7, d'où TES) Sn 


CR | 
On peut identifier le coefficient du potentiel 


par l’étude de l'équation d'état (second coefficient. 


Et 
graphe précédent : ar peur les molécules non 


a Vi + ap? + D pour les molécules 
polaires. Margenau {6] a trouvé un bon accord 
pour un certain nombre de gaz. & 

Dire que l’état d’un gaz peut être représenté. à. 
l’aide de deux constantes caractéristiques, c’est dire. 
qu'il existe une loi des états correspondants. Or on 
sait que cela est près d’être vérifié, mais que ce 
n’est pas rigoureux. La manière la plus simple de. 
rendre compte de ces petits écarts est de prendre 
un ane s variant un peu d’un te à l'autre. 


polaires, et 


4. Énergie de vaporisation. — On peut aussi 
chercher à relier E» à l'énergie de vaporisation. 
Désignons par d la distance moyenne entre molé- 
cules voisines dans la phase liquide, et soit w(r) 
la fonction de distribution du liquide [7] : le nombre 
de centres de molécules compris dans une couche 
sphérique de rayon r et d'épaisseur dr autour 
d’une Do centrale est 


Ar IN” 
4 


w(r)r? dr. = UE 


Aux grandes distances, w(r) = 1. 
L'énergie potentielle de N molécules est 


NE 


as LOTO r2 dr. .. (48) 


C’est l'énergie nécessaire pour séparer les unes 
des autres les molécules du liquide, c’est-à-dire, 
au signe près, l’énergie de vaporisation ; 


Q—RT, Ps - (9) 


en désignant par Q la chaleur latente de vaporisation. 


(20) 


Re 1 Von ne tient pas compte de l’ordre à petite 
_ distance de l'état liquide, on prend w(r) — 0 
_ pour r <rç et w(r) =: pour r>r On trouve 


Line (22) 


éseaux cubique centré et cubique à faces centrés, 


Ndi= FV2: (23) 


: relation (4), ie 
Nrè=V, (24) 
8,=5,60» . (25) 


On peut au Chase schématiser à l'extrême 
la fonction w(r), en:« gelant » la structure réticu- 
_laire du liquide qui n’est en réalité que statistique, 


couche de z voisines immédiates à la distance d, 
puis de couches de 2, z,, ..., molécules aux dis- 
_ tances de . (& est la coordinence moléculaire). 
- Si l’on ne tient compte que des voisines immédiates, 
‘à on obtient, avec d = r», 
Z 
8,— -0». (26) 
2 
Tenons compte des trois premières 


couches. 
Pour un réseau cubique à faces centrées, : 


DSÉCREEN DRE LÉ contribution : 6 Æ,; 
35 — . 6, ds = dy°, Et 0,09€» 5 
gs — 94, d; = dy3, ». OO 
_ Pour un réseau cubique simple, à 
nm, 6, d = Try contribution : 3  Æ,; 
Zo — 12, ds dv», i » 1,39 m; 
z3— 8, .d:= d\W3, » 0,26Æm. 


_ Comme contribution des molécules plus éloignées, 
prenons le terme correctif 


ne p° É Œ SNET Te) 9 
Be f nd En; (27) 


‘avec pour. limite inférieure de l° ie la moyenne 
entre dy ét d, —2d:d,— 2,176; ce terme correctif 


®, 
Les valeurs du rapport 3 sont 


mr 


2 vaut io,01 Ex. 


- done respectivement 


(21) 


Pour des boue oaibattes rappelant les. 


La istanee d est voisine de r», donc, d’après la 


et considérer qu’ une molécule est entourée d’une 


1407 


8,4 et 5,5. On n’obtiendrait 
une relation précise qu'en connaissant bien la fonc- 
tion de distribution w(r). Les données déduites des 
diagrammes X de liquides sont encore trop peu 
nombreuses. Q et w(r)-doivent être pris à la même 
température (en plus de la dilatation, il peut y 
avoir changement de la structure quasi réticulaire). 


Par l'étude de la forme des isothermes avec leur 


modèle, Lennard-Jones et Devonshire [8] ont 
trouvé la relation 
PE Ê 
RE ee 9 
9x 9 ( 7 


La relation (16) serait donc obtenue pour z voisin 
de 11. En utilisant la relation (16) pour toutes les 


0» | 
structures, on trouve une -valeur de pouvant 


varier de 6,7 à 4,4. En fait, ce rapport . souvent 
voisin de 6. C’est une forme de la règle de Pictet 
et Trouton. D’après une règle de Guldberg, le rapport 
du point d’ébullition sous la pression atmosphé- 
rique T4 à la température critique est voisin de 0,6. 
D'autre part, la formule de Cederberg, 


Gr pese 

Q P; P 

DR es po (29) 
Te 


est une bonne approximation. Pour de nombreux 
corps, la pression critique est de l’ordre de 5o atm; 
le point d’ébullition normal est un point corres- 
pondant; la formule précédente. donne alors la règle 


de Pictet et Trouton : Ge 


l’entropie de vaporisation = 
est d'environ 21 cal par degré par molécule; d’où 


> @;, n a} 
7 = 9,9.X 0,6 A (30) 
Remarquons qu'en prenant l'équation d'état 
de Van der Waals, avec les notations Rose 
l'énergie de vaporisation est aussi égale à =-+ Alors 
Te} a = 
UENRP NT co 


Le second coeflicient du viriel, NXTb — a, s’annule 


pour - 


(32) 


2 : ee 
(avec la valeur b — 37 Nri du covolume ). D'où 


O,, D 
ES Tor 


(33) 


valeurs de 7. 75 Si lon s'en 
WE 


tient strictement à l'équation de Van der Waals, 


et l’on trouve pour 


: D fe Tu 
et 5,7, si l’on prend la valeur expérimentale de T. 
te 


e** 


108. 


5. Dynamique des chocs entre M Ce 
Deux méthodes sont employées en théorie ciné- 
tique. La méthode purement statistique des libres 
parcours moyens de Clausius et de Meyer, appliquée 
au modèle sphérique rigide, cherche à suivre le 
mécanisme physique des phénomènes. Ses résultats 
manquent de rigueur, car elle ne fait pas intervenir 
les conditions dynamiques des chocs. Dans l’étude 
de la diffusion, elle doit recourir à l’artifice d’un 
mouvement d'ensemble du gaz, le calcul brut 
indiquant l'établissement d’une différence de pres- 
sion. ( Voir aussi la théorie élémentaire de la diffusion 
thermique que j'ai publiée dans ce journal) [9]. 
L'autre méthode, introduite par Maxwell, analyse 
complètement les conditions des collisions inter- 


moléculaires. Les calculs sont très compliqués, sauf 


dans le cas du potentiel d'interaction de Maxwell 
(fig 1 b avec s — 4). En effet, dans ce cas, la fonc- 
_ tion de distribution des vitesses n'intervient pas. 
_ Dans le cas général, c’est justement le point de 
départ laborieux d'obtenir pour un état non stable 
correspondant à un phénomène de transport la 
fonction de distribution des vitesses Pet correction 
de la fonction de l’état stable. 

Le mouvement relatif de 2 molécules de masses m 
et m, ayant une énergie potentielle d'interaction E(r) 
est le même que celui d’une masse unité autour 
d’un centre de force fixe avec une oi . potentielle 
à la distance r 

Er). nu+m ï 


PREND (34) 


b I No 


p. étant la masse réduite. Prenons des coordonnées 
polaires r, © dans le plan de l'orbite. En désignant 
par p la distance du centre de force à l’asymptote 
de la trajectoire et par g la vitesse relative à l'infini 
avant le choc, on obtient pour le mouvement les 
équations classiques 


RER | 
f w dé DE) 
l +n(s As Er) de 
21 (de 7e AE DEN ven 
Posons £ —$, L'équation de la trajectoire est 
donnée par 


* (36) 


8 2 
RE 
20 
MR ae 

ne m8? 
dr 

pour de —= 0, 
c’est-à-dire pour la plus petite racine f; du radical. 

L’angle des asymptotes est Fe 


6; d el 
ef — 
0 


Le périhélie (ou apside) est obtenu 


— Quand il n'y a 


_(diffusion, conductibilité thermique), la seconde po 


Remarquons que la moyenne de + est un repère 


C2] 


ee Ne 
£ . = 3 : 

IE is [ sin? 9; cos? ®;p dp = | - sin? gp dp. 
“0 à 0 


On les appelle sections de choc, la première p: L 
les phénomènes de transport de masse ou d’énerg 


les phénomènes de transport de re 
Det (viscosité). Dans le cas Rue m 


les différentes sortes de chocs ( — 1,2 — 2,1 — 
Prenons le potentiel d'interaction (5), et posons 


P me 
7e FO ERA 


Alors 


I(5) AN s 

“ap 

1C) 2 7 JP sin 9 

ee “ 29: d(p?) 


de température; en effet [1], 


> Er 


ML —= 


ce D'où, avec (16), 
ÈS 


T=— 


ss. 


J'ai effectué les intégrations graphiques sur les 
formes (41) et (42). Le tableau I donne les résultats 


obtenus pour ni et 15 Les courbes (fig. 3) oi 


log ee ; et log 
même dise) sait une perturbation locale de pente : 
un peu plus forte pour log J au voisinage de + — = 


qui semble supérieure aux erreurs d'intégration 
graphique. 


ité. — [La méthode élémentaire des 
du coefficient de viscosité 


: ne é = 
= Rama, : -(45) 


ee . (A6) 


> 


la Loule nie. 66) donnée plus bas: pour 


= 


sS + 
25 


2 (49) 


suivant < son raisonnement, on trouve pour la molé- 
cule un diamètre efficace 6,7 donné par 


2 
def: A D en iEn (50) 


valeurs obtenues pour Tom (Ag-2) où les deux courbes 


moyens conduit à l'expression 


et 


: Pour de nombreux gaz, la constante de Sutherland S 


Le de molécule, J — ir. Le Donc 


409 


(1) 


Essayons de relier approximativement E%, et le 
minimum Æ, du potentiel rd de Lennard-Jones. 
Pour rapprocher le plus possible les deux formes 
de courbe, prenons T = To et faisons coïncider les 


AE 


sont en trait plein et les parties répulsive et attractive 
du potentiel id en pointillé. Alors, 


bp}; 2E ; F (52) F 
Sp 087... 69 & 


adoptée est voisine de cette valeur. Posons 


+ d'inm 
Sd In ( 


S est généralement déterminée par la valeur expéri- : 
mentale de y à o°C 


3 T I S e 
2 


ST 0 Si 

Quand T varie de o à æ, y varie de r,5 à 0,5; 
pour T=—S, y —1:1. La formule de ‘Sutherland 
représente assez bien les résultats expérimentaux 
et la décroissance observée de y. 

Hassé et Cook [11] ont fait les calculs pour 
un potentiel de Langevin-Keesom avec s°— 4. Fe 
Leur résultat ne peut s'exprimer en termes finis, 2 
mais leur table de valeurs calculées indique une > 
variation de y de 1 à 0,5. Ils ont pris ensuite un ae 
potentiel de Lennard-Jones avec $s = 8 et s' = 4. 
Ils font leurs intégrations analytiquement et 
utilisent des développements limités. Les courbes 
obtenues indiquent aussi une variation de y de : 
à 0,5, mais avec un maximum dans la région 0,4 Te 


10 
ÿ Remarquons qu'avec s = 19 et s = 6 


et qui va s’atténuer encore dans le calcul de n 


(courbe K de la figure 3). D'autre part, je trouve 


cet accident local vers 0,8 T.. Cela-semble en accord 
avec le fait qu'avec une forme de pe plus 


raide (s ou s’ plus grand) le rapport 9 —— augm nte 


Ge. fig. »). 
La formule générale valable dans tous les cas a 
été donnée par Chapman : 


DTA EEE à 
£ A A Es 56 
ASUS OT RE Gi) 
avec Ée RE RTS 
+. I e 
Î= TS (57) 


e est un terme correctif dépendant de s et au plus 
de l’ordre du centième. L'intégrale K est, en posant 


C=jue= sine = (58) 
k=f re (39). 
9 . 
_ Avec l'expression (40) de +, on trouve 
Fe Ô— _. = or (60) 


La courbe de variation de log X avec 0 a même allure 
P. que celles de Z et J. La perturbation de J est étalée 
et à peine sensible. Posons 
ne | 
sav mure A 
c’est-à-dire 
a Fee 
1 = 5 HO. 
Les résultats sont rassemblés dans le tableau II. 
On peut essayer de les représenter par une formule 
empirique. La formule de Sutherland ne donne pas 


une variation de y assez rapide, et en fait, pour la : 


faire cadrer avec les résultats expérimentaux, on 
doit prendre une constante S de plus en plus petite 
au fur et à mesure que la température augmente. 
Cherchons une formule du même type 


(63) 


(64) 


0,4 Te Fa 
(+ _ 


j'ai trouvé Elle n'est pas 


_ pour J un résultat analogue mais moins accentué 


(62) 


ratures ie 0,6- Pa): 
quand T tend Vers. zéro. 


PAL EAU FRE RM 


6 A 70 y 
0,9 0,8 O0 al qu 
0,63 1 FO AS 0,95. 

1 1,0 0,28 78 
1,20 2 0,23 0,73 
2 3,2 0,15 0,6; 

Â 6,4 0,10 06027 


Pour ie molécules bois, En diminue un pet a 
quand la température augmente, à cause du terme 
d'orientation (2). I s “ensuit dans les formule 


précédentes, par. exemple, une diminution de S 
ou S', c’est-à-dire une diminution supplémentair 
de y. La variation de y doit donc être plus rapide. 
que pour les gaz non polaires. C'est ainsi que, 
lorsque la température varie de Te à 2 Ts Y passe 
environ de 1,05 à 0,75 pour l’ammoniac et de 0,96 
à o,73 pour l’éthylène (d’après les résultats de 
TFrautz et de ses collaborateurs) [re]. ; 


7. Corrections Ar ns “+ IE ÿ na de. 
deux sortes. La première est due au principe 
d'exclusion de Pauli. L'erreur introduite résulte 
de l'application de la statistique classique de. 
Maxwell-Boltzmann à un gaz qui suit les statistiques : 
quantiques de Bose-Einstein ou de Fermi-Dirac, 
suivant que la fonction d’onde attachée à la molécule 
est symétrique ou antisymétrique. On trouvera 
des références bibliographiques dans [4] et [13]. 
Cette correction est généralement négligeable; elle 
n'interviendrait qu'aux très basses températures  … 
ou aux hautes pressions. La seconde est due au 


es. hu somme, au Féolenl “Ho 
à la déviation par le choc se superpose l'effet de 
diffraction de l’onde associée à la particule incidente 
par la particule choquée. Cette correction n’est 
sensible que lorsque la longueur d’onde associée 


“YmkT" 


devient comparable au diamètre de la molécule. 
Cette correction n’est importante que pour l’hélium 
et l’hydrogène, mais dans ces deux cas jusqu'aux 
températures ordinaires, puisqu'à o°C les valeurs 


À 
de = = sont encore 0,061 pour l’hélium et 0,069 pour 


et Mohr (41; disons qualitativement que, lorsque 
a température diminue, la section de choc augmente. 
Suivons un raisonnement analogue à celui qui a 
été tenu pour les molécules polaires. Quand la 
température diminue, la distance r, augmente et 
‘la partie attractive du potentiel ne change pas 
(molécules non polaires); cela revient à dire que la 


que y diminue. Cette diminution doit se superposer 
à l'augmentation du cas général. La compensation 
_ est telle qu’on peut réprésenter-les variations de la 
_ viscosité dans un très grand intervalle de tempé- 
_ rature avec une seule valeur de y : 0,65 pour Fhélium 
et 0,68 pour l'hydrogène. Le fait que la valeur de y 
soit plus petite pour l’hélium que pour l'hydrogène 
_ est en accord avec ce que nous avons vu précé- 
 demment, puisque la température critique de 
-_ l’hélium est inférieure à celle de l'hydrogène et que 
les corrections quantiques doivent être du même 


s : À? , EE 
ordre, les valeurs de “ étant voisines. 


8. Diffusion thermique. — Je rappelle l'équation 
de définition des coefficients de diffusion D et de 
” diffusion thermique Dr 


a(Wi— W)= — D grad & + _ grad 7'; (66) 


©» @& sont les titres moléculaires (ou rapports 


: ca LR $ nl = = 
molaires), €, — Se a en D 
} 1 2 ARE é 


_W,, W, les vitesses d'ensemble des gaz 1 et 2; 
W=cW,+aW, la vitesse de convection du 

mélange. 

Supposons W —o; on peut aussi écrire cette 
équation en faisant intervenir les nombres de molé- 
- cules par centimètre cube n au lieu des titres molécu- 
laires c. En désignant par v, le flux des molécules 1 
. par seconde à travers une surface de 1 em? normale 
- à la direction des gradients, on obtient  - 


V4 
HAE Er 


Le = ni 


(67) 


Rise (65). 


D rosène. Sans reprendre les calculs de Massey 


profondeur du trou de potentiel E,, diminue, done 


Cela revient à remplacer Vm par 


Ant 


A avec 
(hi + R)T= const. 
ou LCR Ro) dei = d'u + 7 . (68) 
N=— D gradn; +{[(ni+m:)Dr— nD] Ds 2 (69) 
Reprenant les premières notations, posons 
. = QC Cr. (70) 


(on adopte généralement les appellations suivantes : 


Dr + 
| 5 lappor 
a facteur de diffusion 


Dr coefficient de diffusion thermique, 
de diffusion thermique, 
thermique ). Le bilan du passage des molécules à 


travers un élément, de surface conduit [9], BSE 


- dans le cas €, < C» au résultat 


Vm2 Bad ÿ ma 
GES MEET 


_ A(vm) 71 
Ne (71) 


qui s'écrit, pour des masses assez voisines (cas des 
iS0LopeS); 


2) Va és 


(a, désignant l’expression de a pour des sphères 


rigides élastiques). 
En réalité, il faut faire intervenir les libres parcours 


moyens L, l, différents pour les deux sortes de 


molécules, et au lieu de l'équation (4) de mon pré- 
cédent article, on obtient 


ES des <e; dif 


Ua d 2 Ci 
Uo b Ï C> 


des PE d In T % 


(73) 


Us u> étant les vitesses moléculaires moyennes. 


Va, 


7" Supposant 


les molécules de même diamètre, prenons pour 
expressions des libres parcours moyens [1] 


L = nl 
— mm: 
(nef ) 
HE) 
(TA) 
I 
DR re Re TS M nes NS CS) 
1 Mo 
= (ne V2 + m\/i+ _ 
« m1 
d’où 
l Am < 
none (75) 
1 4m 


On obtient finalement, au lieu de (72), la formule 


3 Am DSTI 7 


li = => 
1 8m 4 Ma+ nu 


(76) 


MINS 


libres parcours moyens, 
rigides élastiques. Pour. 


s'applique à des sphères 
ce modèle, on emploie 


d'Enskog 
105 75 — M; Am 
Up = —> —— = 0,8 , 
IIS Mo + Ma 2m 


<O 


Cette formule est en réalité le premier terme d’un 
Je “ae à M2 — MA 
ne Dent en série de puissances de —— 
Mo + Ma 
voisines. Chapman avait aussi donné une formule 
se réduisant pour © < € à 


17 Mo— ma Am 


do = — 0,62 EL 


27,45: \Mo + PU A 


(78) 


Posons pour l'expression de a, dans le cas d’un 
modèle quelconque, ; 

M2 — MA 
Mo + Mi 


& = As RT= 


(79) 


Avec le modèle 1b de potentiel répulsif en r°*, 
les calculs de Chapman et d’'Enskog donnent 


Rr=— s 


(80) 
où C(s) varie de 0,9 pour $ —12 à 0,8 pour $ — 4. 


_ pour des molécules maxwelliennes (s — 4). 

Clark Jones [16] a étudié différentes formes de 
_ potentiel rendant mieux compte des interactions 
_ moléculaires : potentiel de Sutherland, potentiels 
_ de Lennard-Jones avec s quelconque, s' — 2 (valeur 
_ qui entraîne une simplification des calculs) et 
avec $ — 8, s — 4. En attendant de reprendre le 
calcul complet pour s — 12 et s' — 6, il a choisi 
ces formes particulières pour lesquelles de longs 
. calculs préliminaires avaient déjà été effectués 
en vue de l'application à la viscosité. Avec ces 
différents potentiels, R7 dépend de la température, 
Dans le premier et le dernier cas, Rr s’annule pour 


c’est-à-dire, en utilisant la relation (16) obtenue 
avec s — 12 et s' — 6, pour la température 


T'ELETe (82) 


Les valeurs positives de R7 obtenues sont ün peu 
trop faibles. L’allure de la variation de R7 doit se 
retrouver dans toutes les formes de potentiel de 
Lennard-Jones. 


9. Diffusion thermique dans le cas de Lorentz. 
HE — C'est le cas d’un mélange de deux sortes de 
molécules r et 2 où les molécules r sont rares et 
beaucoup plus légères que les molécules 2 (dans les 


Cette théorie élémentaire, utilisant la méthode des : 
généralement la formule suivante tirée des calculs . 


cule 1. 


ét n’est donc valable que pour des masses assez 


“OA pose 
ctm=p=m 


Le phénomène de diffusion thermique s’annulerait 


_ Faisons d’abord le calcul dans le cas d'un . I , = 


(81). 


La fonction de a des vitesses 
à 


où fo est da innetion pour l’é état stable 


e—jmnu? 


Pour les phénomènes de transport de masse 
d'énergie, l'expression de d est, en prenant con 
axe des x la direction commune des gradients, 


pour simplifier l'écriture, 


une surface de : cm? est 


re z) erimui É < à IAE RS 
dr 5 3 | 


ue 
LU — — 
Ë ul | ri de 


1e dus du. dus — 


+ 


répulsif en r*. Due ce Cas, 
= S—k 


Toul — na AUSE 


En identifiant les expressions (69) et (0) de vw 
on obtient 


D = —— s 
3Vx : È 
: na A(ÿm) 
= F î 
Sr(2+ 2)r(s+ 2) 
2 $ S 
d da = : : 
F (2+ à 
mais 
r(s+2) = b+2)r(ee 2). 
S > ; 
d’où 2 
d = 5 Ë 
j à 2:S = 
c’est-à-dire 
du = =, D, 


is da, est plus petit que = valeurs 
nt les formules (76) ou (77) appliquées 
em au cas M Æ Mo ae des rapports 
de Ar et ie 


4 (jm)  … du - 
j eo — ————— 99 
2: 3Vx 7 ë I re (92) 
= É e—imé y du 
DS 4 0 I + 
Re ao (93) 


eines du 
Es FE 


Ti 
… 267 


: est le même pour les divers 
d'interaction Pi 2 ie, ONE peut 


EE 
TT + En réalité, ie serait une moyenne 


< les un critiques des deux gaz. Le 
facteur de diffusion thermique «a; s’annulerait done 
2 voisin de 0,55, ( SA ure Dour une tempé- 


PASS 
c sx 


Pi ST. 


ette_ 1 est sensiblement égal 


que as nie à peu près pour la même ne 
re 0,9 T., mais que le coefficient b de la for- 
mule (79) soit supérieur à la valeur a; du cas de 
HAN: Dans le RAT III sont portées les valeurs 


10. Approximations et cas des gaz de masses 
moléculaires égales. — On peut chercher dans 


onnement de la théorie des libres parcours 
moyens. Il faut alors introduire un diamètre molé- 
culaire variable avec la température. Pour tenir 
compte de la variation des libres parcours moyens 
avec la température, remplaçons l'équation du 
bilan de la diffusion thermique (73) ou 


= (die - din r) (95) 


formule à peu près équivalente à la relation 


(à Dee lan dinT — à Ins:) 


AR qu 
. 


en supposant, comme précédemment, les diamètres 
des molécules r et 2 égaux. Cela nous fournit, dans 
le cas qu Re répulsif en r*, la valeur approxi- 


(ame? dnT—dinst), (96) 


mative © fie Rr. Dans le cas général, ce ne sont 


pas en toute rigueur les mêmes sections qui inter- 


viennent dans les phénomènes de transport de 
masse et d'énergie et dans ceux de moment; mais 
étant donné l’analogie d’allure des courbes I, J, K, 
on peut prendre pour dInc? le résultat donné par 


la viscosité. Il faut donc, dans l'équation (95), 


remplacer : par = — 0, c’est-à-dire multiplier le 
résultat (76) par 


1— 20 = 2(1 — y). ee (97) 
On obtient 7 
= (NET (98) 


LP 


Rr=1,710— 7%) (99) 


obtenue par Clark Jones dans le cas du potentiel e 


répulsif nr Ee les formules (77) et (79)]. 
Les valeurs de en) sont portées dans le 


tableau III en regard des valeurs expérimentales 

de b obtenues par’ Nier [17] dans les séparations 

isotopiques du méthane et du néon. La tempé- 

rature T indiquée est la température moyenne 
de Brown [18] 

PLAT B Ts 

Po EE 

où T4 et Ty sont les températures extrêmes. Elle 

ne diffère pas beaucoup de la moyenne géométrique. 


TaBLEAU II, 


129 
Nes à ni A500 D o 
408 0,63 0,94 0,47 
(extrap.) 


Dans la diffusion thermique de l’ammoniac, 
j'ai trouvé [19] un effet isotopique beaucoup plus 
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TT PARLER ï 
AO A Re 


RÉ es à ù Fa +. : Prec a ah \ 
LIT EC rather 


“ENS faible que pour le méthane dans les. mêmes condi- 
5 D eo tions (T4 300, T5 = 550), mais positif : les 
RU 0 : molécules lourdes se concentrent bien vers les 


Fig. 42 


régions froides. Le potentiel répulsif en r—* conduit 
_ pour l’ammoniac à s — 3,6 et fait donc prévoir 
un effet négatif. Ce potentiel ne convient donc pas 
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action moléculaire que la diffusion 
L'accord qualitatif obtenu avec le potentie 
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de Wall et Holey relatifs à la diffusion thermique 
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figure 4 présente. À masse égale, les molécules 
qui vont vers la paroi froide sont celles qui ont. 
plus grand diamètre (& ou r,) et aussi la plus grande 
polarisabilité, donc le potentiel attractif le plus 
grand, les potentiels d’ionisation ne jouant pas le 
rôle prépondérant et étant du même ordre de: 
grandeur. D’après la forme des courbes, ce sont. 
donc aussi les molécules qui ont le trou de potentiel 
d'interaction le plus profond, c’est-à-dire la tempé- 
rature critique la plus haute. Cette conclusion 
est vérifiée pour les paires de gaz suivantes : 
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. de coïncidence. & 


Nous nous proposons ro ce travail 1 étude de 
has Jhoations pouvant s’exercer 


ue neutrons 


: u ES botentiel de Coulomb, tandis que les inter- 
actions s’exerçant entre protons ou neutrons par 
l'intermédiaire du méson conduisent à la théorie 
je Les RUSISAIres Nous allons montrer ici que 


nt d'interactions permet de représenter les 
interactions massiques correspondant au potentiel 
de Newton. 


“1, “La particule de spin - . — Les particules 


DL et ; RTS 
de Spin, >> 53” protons, neutrons (nucléons}ou électrons 


seront cents en l’absence d'interactions par 
les ondes d,(X, #) (i = 1, 2, 3, 4) solutions d’une 
équation de Dirac 


h 1 d, RNA 
nous rar sp mac | 
| Les As %4 Fe, quatre matrices de Dirac telles que 


(a J = 1, 2, 3, 4). 


CR ajay oja= 2j; 


_ générale de 1’ interaction liant les particules de spin 


L 
actions s ’exerçcant par l'intermédiaire des potentiels Tlongitudinaux de la particule de spin 2 . 


ee | SUR LES INTERACTIONS ENTRE PARTICULES MATÉRIELLES 
.  SEXERGANT PAR L'INTERMÉDIAIRE DE LA PARTICULE DE SPIN es 
Se Par Gérarn PETIAU. 
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* Sommaire. > Par un raisonnement analogue à celui par lequel on retrouve le potentiel de Coulomb dans 


RS Cr LE e 3 De é 
…_ la partie statique des interactions entre corpuscules de spin ;— et photons, nous montrons dans ce travail 


4% 


comment la théorie de la particule de spin 2 22 permet de retrouver en considérant cette re 


comme Dent d'échange d'énergie et de HUAMULE de mouvement entre des particules de spin — les 


interactions massiques statiques correspondant au potentiel newtonien. Après avoir établi une Lois 


k h 
— et 2 — le calcul est particularisé au cas des 


nd 2T 


et l’on Ë 


ua 


montre que l’élément de matrice d'interaction calculé dans le cas statique est équivalent à l élément 
correspondant du potentiel newtonien corrigé “DAe un terme exponentiel et complété par deux termes 


7 


Nous adoptons ici la représentation usuelle 
= — ps; (6p6g=— 107; PpPya=— tr) 


les co, ou p, sont les 6 matrices. de Pauli. 
Nous introduirons également le tenseur métrique Ju 
tel que es 


Ap = Sp; Cr 


800 —]; Spp= 1; 


et nous poserons 


UNE le Up =— UP = àp; O = % ; 
Uop = — UP = 1); Ung = URT = Ta u 
Posant Fe 
) at ) d 
DOVE LENS OR= 
e.ot: CEE 


l'équation (1) s'écrit 


h. | 
. JAutThi= mocoŸ; 
DTA 


D Pi, P5:, P:50), (2) 


RESTE 
Se 0 ! CRE NE: Pi iris 
XL C { 5 Ê + 1rtolasts 


+pocal at [Prin (3) 


Ë Co 0 PONT RE 
Die mr dp= = dde ds, = 1,234): 
( t c dt? D 0€ 9 55 V4 RS D) 


D, 


_ du second ordre 


On peut montrer que les solutions de cette équa- 


tion se répartissent en deux classes. 


La première classe comprend les solutions de (3) 


qui sont également solutions de l'équation 


/ 
4% D SRE ART 
CP; à ER WE H0 C* Din 


On peut montrer que ces fonctions constituent 
l’ensemble des solutions du système des quatre 
équations de Dirac simultanées 


h - 
ar dy (ab Mme De, Mb = HO c( 0) lim D; RATIO (4) 


27RL 
(r=1, 2; 3, D 


Ces solutions ont été étudiées en détail dans le 
livre de M. Louis de Broglie « Théorie générale des 
particules à spin » (Chap. XI, p. 182). Elles se répar- 
tissent suivant les états du spin total et caractérisent 


un état de spin total », trois états de spin total : 


et deux états de spin total 0. En particulier on peut 
montrer que les solutions de spin total 2 sont les 
fonctions ©; ; 
rapport aux indices dy, ds, les de 


Les solutions de la seconde classe qui sont solu- 
tions de (3), mais non de (4), satisfont à l'équation 


TL RO De 
h2 HE Doi iast». 


Pia = — 


x 


Ces solutions correspondent à des états de spin 


h 
total = OÙ ZéTo, d’un HR de masse 2 
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La constante Bo masse propre du corpuscule 


ï 


représenté par l'ensemble des solutions de (3) 


__ doit être vraisemblablement sinon nulle, du moins 
très petite et du même ordre que celle du photon 


si l’on considère que l'équation (3) réalise par l’en- 
semble de ses solutions une théorie unitaire des 
champs gravifiques et photoniques. 


3. Ondes planes de la particule de spin 2 - 


nm 


— Si l’on cherche une solution de l’équation (3) 
dans le cas du spin total maximum, sous la forme 
d’une onde plane 


se: > 
RE ALES ñ —ikæ 
Pit —dlithhe ? 


on trouve pour déterminer a, l'équation 


) lastgsda 


avec 


complètement symétriques par. 


Rap R) = as Ra RE af; 


ee a. 
X [Æ — hÙ) ] [Æ — Al) ln. ska 
X Amimimsm ele) | de 
+  [k+ AMI(E + A@)] 
x [k + ABLE + Aie rs mn 
: ; > - 
Æ Br ….,m, eitte+kt)s 


ll 
Pitiiks ki = 


Àest une constante de normalisation. Les 
Bnsmsms sont deux ArSUpes de 16 € 
positive, Les FR sont Sous nn ce que n 
supposerons par la suite et il reste 


Dirk, = MES pu) Mk Ro] — 2 
T4 — RE DEA — A0) a 


X At mm eRR %), 


At ee 

P* FE ( Dors . ar 

d’où 
9 RE) _. 4 x Le BP De 

162 ÆX 0 ( k+ ko )' En,memam, ( AA }nssmasmsn © = 2 


facilement que certaines combinaisons hits de 
des D; 11, Se comportent comme des tenseur 
Soit @’ la matrice telle que : 


(nous désignons par A* la matrice telle que = Le #7 
avec la représentation des « choisis, on voit que 


RL TO 5 A —= — 09 p2. 
L'introduction de @’ donne 


(R'u? DZ == R'ur, 
(R'uPIYE= R'urT 


CRUE Ru, 
(R'ur)ÿ+= Ru, 


et, par suite, ces matrices sont symétriques. 
Avec l’aide de ces matrices, nous pouvons « 


& Ru), (@ ut de 
pare Dre, 
,, (R' ua Co ur jai, D ae 

DRE era ab Ar 


blostgsts 


= 


“consta ate de proportionnalité qu ’avec 


n 0, Pas, = Dot 
des: Fonte — Par rapport 2e 
ction de propagation, les ondes tensorielles 
écédentes peuvent se répartir en ondes longitu- 

s L -transversales et transversales. 
direction de. propagation étant 
2£. les Co ce se 


0,0 . ni De2,0,. 03,3; 03,05, 


On. voit alors que dans le « cas d'une onde “ha 
ndeurs ne on que de la combinaison 


ae = Au == oe = Are A ESS 


I 


AU RNTAI = Re 


=. & ne . _ 
É VAS ko le 
qu 2C#l#|P (11) 
D Voie - 
— 
$ VAS ke | 


P = eiltkctk2) 


(9) 


Ë no 1lN [al CAR 


; - 1 Ÿ Î 7 l 
Eee — ï > D On 
Pe As ass 

(10) ee + | 
- H,, est l’opérateur hamiltonien de la 


2 & T | 3 
h-\-d 1 ; le 
JE —°C me CANAL 4 De. Caro & 
© . a 
Er SRT) OL 4 >. RER De none HSE 
: P Le Û 
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tudinaux : 
RE ko 
P03,3 — 2CKkoiP 
Dos 0 NE : 
ES 5 VAi ke 


6. Le système corpuscule de spin —— Corpus- 
É ( É 


Pnee Se = Fe 
cule de spin 2: ee — Nous représenterons le système 
44 ’ k- 


total en interaction par l'équation d’ondes 


2 


AE PP ES 
DRE Fe A btrtmiasts 


Hpn+ SC Hi bilatbse (12) 


= Here 
particuie 


de spin a # celui du graviton LEZ 


h De 
Hp= cc 3 0, Xp + mes}, 


H; est un à opérateur d'interaction que nous écrirons 
Cote 80% B+ mtap,y 02BY 
S 5 > | < 
7 rys(R —7). (43) 


Les 1,8, {184 tlas,yô sont des opérateurs de la 


; ee & ARE 
théorie de la particule de spin ne 


Nous poserons 


. RE - MC 
0 : Las © 


mo C 


= bp —, sago = tag mL (14) 


1,8 est le tenseur de Tétrode et l’on voit facilement 
que l’on a 


la 1 0, 


» 


I RENTREE à <e 
lub,y = à (=) (ca Hy0Br = 0Btyde) 
h 
Fan (td aus) | 


h no 
DU RAA EU de w + w0 
(=) 6 | Sox (28 8 “ 


— gpy (220 + 009) ] (45) 
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et, par suite, 


las? = 0; 
1 h \? 
pre = dy llyg08— dGU;5 da 
oi Dr \2Tt De ENT . 
+ djuagds— 05 Uasoy) 
“ < STE —> 
M € 


—g (. Barth #py60 
On a alors par contraction 


270 C 


Cap, v8 = = dust. 


Les b/ sont des PRESSE HA de la théorie 


de la particule de spin 2- 


2 Pi 
( baby da Mas ee LNy 
I 
Fe l (utuÿ+ u? us) o,—(utut + ut ee 


HAUT RUE : A) EG (ue uR+ u}ut)oiws], (17) 


lié à par, Ja relation 
PA > z 
7. à mMocn = l. 


7. Étude des éléments de la matrice H (!). — 


- En. l'absence d'interaction, les fonctions propres de 


_ l'hamiltonien du système corpuscule Le spin . 


ch 
_ corpuscule de. Sp 2 . sont repré par des 


produits WP} 1,1, d'une onde plane du corpuscule 


de Dirac par une onde plane du graviton. Si le 


_ corpuscule de Dirac n’est pas Pie mais dans un. 


- état stationnaire d'énergie En= * ee he les Fi sont 


ee de la forme 


gr = La, un (A, Y A) eiKmct 


és ee ee ee Rx) dk . 


superposition d'ondes planes. 

Dans cette étude nous considérerons. les particules 
de Dirac non engagées dans des liaisons et nous consi- 
dérerons les transitions du système d’un état de 
mouvement y” = ddl}, produit d'ondes planes, 
à un état dans lequel le graviton est annihilé, c’est- 
à-dire représenté par la fonction constante 


0 2 
D; i = (4 RT! Jus ( ne DTA ë, 


et, par suite, 


(ES 
Hyse 


ii (R'a; Vic ( R'a ie 


Dans cette transition que nous représenterons 
par le symbole (0, m'; l, m), l'élément de matrice 


HE Lin) fl rl R)a der \ 


+ gBèlay— SadlBy). (16) 


l'et m sont deux constantes indépendantes, n est 
l 


tion (18). se ee à 


Or, on voit Hplenieus que Le 
; CH*(b«b)D = Dab, 


\ 
: 72 


D'autre part, nous poserons A 
Ne sde > 
qe )x #) (tab hi; 2 ns TL mn) PR JES 


Nous aurons alors 


Ho: ;Ln) ele  - Va ( 


al 
+ nTas Cr )eg, Cr). ar. 


Nous calculerons explicitement le terme 
dans le:-cas- de es par. les cn 


8. Calcul plots de Hi. Fe. Ÿ: 
cas des ondes longitudinales le terme en Lde ré équ: 


He im = 0 To, Das ©: 
oo 


Tag= Ta 


Posant 


et introduisons pour les l'E les potentiels long 
tudinaux normés (11), il vient 


QUES = Lie a | rs, M k ps Ty ll nl : 


Vki ko 
mue ete | al Jetta ea 


(Ke Pr m) = yln' a a D 


, 
F > 


tels que 5e 
MA EE (Ba B ) UN) cn kan Ke) L 
>t è À ! l x 


nous obtenons 


JOLI EER HIS 
HO HR EP errE== £ 0 | k|+%, :k|k 
Vkÿ ko ) f 0,5 1" 
0 b : 
; "# mc. : 
+ dk? + — io | 


(21) | 


>< ke RER 071 Fr : 


eo or. k = OS 


C dm lorsqu’ il ya conservation de la quantité 
de mouvement. Mais alors l’ exponentielle dépendant 
du temps n'étant pas nulle, il n’y a pas conservation 
de. l'énergie. Il faudra donc considérer au lieu d’une 
transition directe, deux transitions successives en 
passant par un état intermédiaire et shine la 
formule de PArauen fe 


HVH) L 

Hip Hp) PT, 3 
_ Dans un premier processus, la à particule (1) passe 
‘de l'état (01) (Ko, A) à l'état (1)(K3, K3) en cédant 
1. + 

Ja quantité de mouvement k à l’onde longitudinale 


Ds 
du graviton qui est émis avec (ke, k). Celui-ci est 
ensuite absorbé par la particule (2) qui passe de 
l'état (02) RE) à l'état 10 _ 4e 
on a alors 


SRE 
re K; = K: ut EE = ka. 


te. processus global conserve l'énergie, c'est-à-dire 
que 
Ki + Ko2 = Ki + K> ou 


? 


pe — Ki = K;3 — Kyo = K 


es qui définit K. 
_ D'autre part la relation de la théorie de Dirac 


(nim) © 


der = 0 
do 108 + 03 73,8 


= 0, 


En Km 18 1e en )erEE 0 


a > 
HR CA Kn— Km =—K%, 
KB,8= KG, B. 
_ La formule (21) devient alors 
É 2C1 
 Hf0:1,0 = Mars Re ; 
Le La LE PEL En 
>: 
TE Se eilku—Ki+hct, (23) 


à _ Pour le processus global, (23) reporté dans (22) 
nous donne en remarquant que 
hc 
EEÈR= ae ŒKi— Kp) 
K;= Ko1 + Ko, Kp = K1+ Kk + Ko, 


hoc 
Ei— Ep= ER), 


4149 


X (@9,6)1, + — 


3 Hiouu | 
cf(lexe ee) 
Fa 


Mo C Fe I 
: ire (29) 


X (D ,0 )2,02 + 


Le même échange global d'énergie et de 
quantité de mouvement peut avoir également 
lieu d’une façon inverse, la particule (2) passant 

>: ne = + ee 
de l’état (02) (Re: K.) à l’état (2) Fe k&.) en: 
cédant la quantité de mouvement — k au graviton 
qui est ensuite absorbé par la particule (1) qui passe 


de l’état (Gr) LE Ki à l’état (K, Ro d'où 


Ÿ 


> Sie re 
RU Koh Ko1 — = K1, 


Koi — Ki = FR ee F 
Ki—K» = (Ko1 + K52) ET Ho ke + Ko) =— (K re k). 


On a alors pour »H) une expression analogue 


à celle de a) dans laquelle le facteur — 7 
à 1 SE 
est remplacé par res 
L” élément de matrice global est alors * 
HO = G)H 0) + (0H) 
16C2l?r I 


ne Kihe. (K?—%), 


*< FERRÉ KK + Gr) 
pl 


X (Bo0,0 1,01 + TE (o ha, «| | 
> À 
«tn 0 
A If 


X (D ,0 )2, 


&(e | (923) 


Or en première approximation, on peut négliger K 
devant k, ce qui revient à négliger la perturbation 


2 : h 
du mouvement des particules de spin ; par les 


gravitons, et d’autre part si leurs mouvements sont 
suffisamment voisins du repos, on aura 


He ' HR 
Goo — 57 Loi + Ai] Go R se 2 Ko(1)1,01 © 2 74 C(1)1,01- 


De même, on pourra poser, 


& (1h, 01: 


(CE 01 © 


D'autre part, avec la valeur de C? — . > il vient 
() 


Gr > 
Ajhe Xi 
et il nous reste ; 
20 + 
(Hope 
: Kÿp 
*< Êr k le + fm MaC?(1)4,01(1)2,02 (26) 
Or, 
- : > |2 kè 2 & [2 ; 9k$ 
d’où 
18/2 Mmiy/Mac? 
HUE Re Giouis 
4 l2 ga k 
Or on sait (voir par exemple L. de Broglie, Une 


nouvelle théorie de la lumière, t. 2) que le terme 


; AU) o1 (202 est l'élément de matrice corres- 


5 
e—ko lie L. 


pondant au 1 potentiel 7 


ot 


correspond. au potentiel Ô (| rl) 


à A 212 Û I 
De même l'élément de matrice Chow (he. 


donc le même qu 
Re potentiel ë 


(= 


de Do ne peut à Un aux he gra 
tances et par ailleurs see considéré dans les 
cosmologiques. Et 

ee _potentiel est ‘complété par de 


dans les théories quantiques + | 
La constante l est déterminée à 


par la formule 


87% 8x 


Pres 


9 © 9 


Pour le terme en “a Gi) (ho. Un raison- taires. Pa ailleurs cette théorie Aa ul Ée 
‘sentation quantique du formalisme de la th 
nement, analogue à eh fait par M. Louis de Broglie  }; relativité générale en abandonnant postulats 
nous montre qu'il correspond au potentiel — a"(fr. DA ), de celle-ci sur la structure de ‘espace et ens L te 
ltoncon à ( | : ) étant lat aux postulats de la théorie de la relativité res re 
a fonction telle que Il semble donc que cette théorie soit appelée : 
“ 2 (A LE ar = f" jouer un rôle très important dans la représentation 
+ sa 4 er (o). générale de la structure se des interactior 1 
rielles. * : as _ Ma 
ï rs .. Manuscrit reçu le 19 mars 194: 
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